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ABSTRACT 

We present an analysis of photometric and spectroscopic observations of the symbiotic 
binary system RX Pup with the aims of developing a reliable binary model for the 
system and identifying mechanisms responsible for its spectacular activity. The binary 
is composed of a long-period Mira variable surrounded by a thick dust shell and a hot 
^ 0.8Mq white dwarf companion. The hot component produces practically all activity 
observed in the UV, optical and radio range, while variable obscuration of the Mira by 
circumstellar dust is responsible for long-term changes in the near-infrared magnitudes. 
The observations show RX Pup underwent a nova-like eruption during the last three 
decades. The hot component contracted in radius at roughly constant luminosity from 
1975 to 1986, and was the source of a strong stellar wind which prevented it from 
accreting material lost in the Mira wind. Around 1988/9 the hot component turned 
over in the HR diagram and by 1991 its luminosity had faded by a factor of ~ 30 with 
respect to the maximum plateau value and the hot wind had practically ceased. By 
1995 the nova remnant started to accrete material from the Mira wind, as indicated 
by a general increase in intensity of the optical continuum and H i emission. The 
quiescent spectrum resembles the quiescent spectra of symbiotic recurrent novae, and 
its intensity indicates the hot component must accrete as much as ~ 1 per cent of the 
Mira wind, which is more or less the amount predicted by Bondi-Hoyle theory. The 
earliest observational records from the 1890s suggest that another nova-like eruption 
of RX Pup occurred around 1894. 

Key words: stars: individual (RX Puppis) - stars: binaries: symbiotic - stars: novae 
- stars: mass-loss - circumstellar matter. 



1 INTRODUCTION 

Despite numerous papers describing the multi-frequency be- 
haviour of the symbiotic binary RX Pup (HD 69190), and 
several attempts to determine physical mechanisms respon- 
sible for its spectral peregrinations, it remains one of the 
most puzzling symbiotic systems ever found. 

Based on the peculiar emission-line spectrum, as well 
as irregular photometric variability, RX Pup was initially 
placed among the irregular variable stars related to r] Car 
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or R CrB (Pickering 1897, 1914; Payne- Gaposchkin & 
Gaposchkin 1938). In 1941 the star was characterised by a 
continuous spectrum confined to the yellow and red region 
and very intense high excitation emission lines of [FeVII], 
[NeV], [FeVI] and [CaVII], similar to those observed in 
CI Cyg and AX Per, which led Swings & Struve (1941) 
to classify RX Pup as a symbiotic star. The presence of 
an evolved red giant in RX Pup was not established until 
the 1970s, when M-type absorption features were observed 
for the first time (Barton, Phillips & Allen 1979, Andrillat 
1982), while extensive photometric observations in the near- 
infrared demonstrated that the system contains a long pe- 
riod variable with P ~ SSO"* (Feast, Robertson & Catchpole 
1977; Whitelock et al. 1983a, hereafter W83). Direct evi- 
dence for the binary nature of RX Pup has been provided 
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by the detection of a hot UV continuum source with the 
International Ultraviolet Explorer (lUE) satelhte (Kafatos, 
Michahtsianos & Feibelman 1982; Kafatos, MichaUtsianos & 
Fahey 1985, hereafter K82 and K85, respectively). 

Because of the large volume of observational data avail- 
able and the relatively recent onset of activity, RX Pup has 
been regarded as a key object in unravelling the nature of 
symbiotic Miras; especially as it is one of the brightest sym- 
biotics at radio and IR wavelengths. The models proposed 
to provide the hot component luminosity and to explain 
the observed activity involve - as usual for symbiotic bi- 
naries - either thermonuclear runaway on the surface of a 
wind-accreting white dwarf (Allen & Wright 1988, hereafter 
AW88) or variable disk-accretion onto a compact star (Ivi- 
son & Seaquist 1994, hereafter IS94; see also K85). 

In the AW88 model a low-mass white dwarf slowly ac- 
cretes hydrogen-rich material from the wind of a distant 
Mira companion until a thermonuclear runaway occurs. Such 
a flash, especially for a low-mass white dwarf, can last for 
many decades; the hot component of RX Pup now has a 
luminosity close to the Eddington value. The observed peri- 
ods of low and high excitation were postulated to be caused 
by variable mass loss from the white dwarf, while the low 
reddening reported during the low-excitation stage was at- 
tributed to the close proximity of a shock front (where the 
hot component's wind collides with the Mira's wind) to the 
Mira. 

IS94 proposed an alternative scenario in which, during 
some portion of the orbital period, the Mira fills its tidal lobe 
and is gravitationally stripped of its dusty envelope; as a re- 
sult the circumstellar reddening is significantly reduced. An 
accretion disk or torus forms causing the UV radiation of the 
hot companion to be reprocessed. As in the AW88 model, 
the bolometric luminosity of the compact star remains un- 
changed as it goes from high- to low-excitation phases. In 
this model the fiattening of the radio spectrum and the dis- 
appearance of the high-excitation emission lines result from 
cooling of the radiation field in the orbital plane. 

Unfortunately, both models have some significant diffi- 
culties when confronted with existing observations. In par- 
ticular, both AW88 and IS94 assumed that RX Pup drifts 
between two extreme - low and high-excitation - states 
without major changes in its bolometric luminosity. In 
fact there are major differences between the low-excitation 
phases reported in 1949/51 (Henize 1976) and the late 1960s, 
when the optical continuum was very faint (mv ~ 12.5, 
mpg = 13.5 in Feb 1968; Sanduleak & Stephenson 1973) 
and the low-excitation phase observed in 1972-1978. The 
latter was characterised by a significant increase in the op- 
tical brightness (V'~9 — 8, 5^9 — 10) which strongly 
suggests intrinsic variability of the hot source. 

This paper contains an analysis of photometric and 
spectroscopic observations of RX Pup made with the ob- 
jective of developing a reliable binary model and identifying 
the causes of its unusual variability. A description of the 
database is given in Section 2, it is analysed in Section 3 
and the results discussed in Section 4, while Section 5 con- 
tains a brief summary and conclusion. 



2 OBSERVATIONS 



Table 1. Journal of observations 



Date 


JD 24... 


SWP 


LWR/P 




1978, 


Sep 18 


43770 


2684 L 


2395 L 




1979, 


Jun 30 


44055 


5680 L 


4922 L 


9.2 


1979, 


Jul 17 


44072 


5835 L 


5081 L 


9.3 


1979, 


Jul 19 


44074 




5106 L 


9.3 


1980, 


Apr 15 


44345 


8762 L 


7505 L 


9.7 


1980, 


Sep 21 


44503 


10189 L 


8856 L 


9.8 








10190 L 












10191 H 






1981, 


Jun 11 


44767 


14239 L 


10830 L 


10.2 








14240 H 


10832 L 




1982, 


Mar 22 


45051 


16597 H 


12837 H 


10.4 








16598 L 


12836 L 




1983, 


Oct 30 


45638 


21402 H 


2177 H 


11.0 








21403 L 


2176 L 




1984, 


Mar 11 


45771 


22461 H 




11.0 








22462 L 


2924 L 




1986, 


May 8 


46559 


28270 H 




11.2 








28271 L 


8160 L 












8161 L 












8162 L 












8163 L 




1986, 


May 9 


46560 


28278 H 


8170 H 


11.2 








28279 L 


8164 L 










28284 L 


8165 L 












8166 L 












8167 L 












8168 L 












8169 L 




1987, 


Jul 3 


46980 


31285 H 




11.3 


1988, 


Oct 24 


47459 


34594 L 


14309 L 


11.5 








34595 H 


14310 H 




1989, 


May 20 


47667 


36307 L 


15552 L 


11.7 








36308 L 






1991, 


Mar 29 


48345 


41209 L 


20008 L 


12.2 


1993, 


May 27 


49135 


47757 L 


25611 L 




1993, 


May 28 


49136 




25615 L 





2.1 Ultraviolet spectroscopy 

Various observers acquired ultraviolet spectra of RX Pup 
throughout 1978-1993 with lUE. These spectra covered a 
significant part of the last phase of activity. Table 1 sum- 
marises the relevant parameters for these observations, in- 
cluding the Julian Date (JD), the spectrum identifier {L 
and H indicate low and high resolution data for short wave- 
length - SWP - and long wavelength - LWR and LWP - 
exposures) and visual magnitudes derived from FES counts. 
Most of the lUE data, in particular the low resolution SWP 
spectra taken prior to September 1980 and after May 1986 
as well as practically all LWR/LWP spectra, have not been 
published. 

The low resolution lUE spectra were transformed into 
absolute fiuxes using standard calibrations (Holm et al. 
1982) and the LWR spectra were corrected for sensitivity 
degradation (Clavel, Gilmozzi & Prieto 1986). The high res- 
olution spectra were calibrated as suggested by Cassatella, 
Ponz & Selvelli (1982). Table 2 lists emission line fluxes 
derived by fitting Gaussian profiles as a function of Julian 
Date; these estimates have errors of ±10% for strong emis- 
sion lines and ±25% for weaker lines. Examples of ultraviolet 
spectra are shown in Fig. 1. 
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Table 2. Ultraviolet emission line fluxes in RX Pup in units of 10 erg cm ^ s 



JD 24... 


44055 


44072 


44345 


44503 


44767 


45051 


45638 


45771 


46559 


47459 


47667 


Nv 1240 










0.3 


0.5 


0.7 


0.9 


1.6 


0.6 


0.7 


O I 1304 


0.8 


1.2: 


1.0 


0.7 


0.8 


1.3: 


0.3 


0.3 


0.3 


0.1 




Si IV 1394 


0.6: 


0.5: 


0.7 


0.7 


0.8 


0.8 


0.7 


0.7 


0.7 


0.4 


0.3 


Oiv] 1403 


0.5: 


0.8 


1.0 


0.8 


1.6 


2.6 


2.2 


3.5 


2.2 


1.1 


1.2 


Niv] 1487 


0.9 


1.4 


2.7 


2.4 


3.4 


3.5 


3.1 


3.3 


3.1 


1.3 


1.3 


C IV 1550 


7.2 


8.5 


10.7 


10.6 


11.7 


11.3 


9.5 


8.7 


8.5 


6.5 


6.2 


[Nev] 1575 






0.3 


0.3 


0.3 


0.3 


0.4 




0.3 


0.3 


0.2 


[Neiv] 1602 






0.6 




0.4 


0.6 


0.5 






0.5 


0.4 


Hen 1640 


0.7 


1.1 


3.5 


3.9 


5.3 


5.6 


6.0 


5.8 


5.1 


4.1 


3.5 


Oiii] 1664 


2.3 


2.9 


5.9 


5.7 


5.6 


5.1 


4.0 


3.8 


3.1 


1.5 


1.1 


Niv 1719 










0.2 


0.2 




0.3 


0.2 


0.1 




Niii] 1750 


1.8 


2.7 


3.0 


3.6 


3.7 


3.4 


2.7 


2.5 


2.0 


1.2 


1.1 


Sin 1818 


0.2 


0.3 


0.5 


0.5 


0.4 


0.5 


0.5 






0.2 


0.2 


Sim] 1892 


2.6 


2.9 


3.7 


3.7 


3.8 


3.6 


2.6 


2.4 


2.0 


1.0 


0.8 


Cm] 1909 


4.1 


4.5 


8.0 


8.5 


10.0 


10.1 


8.3 


8.5 


5.9 


3.8 


3.4 


He II 2511, Foil 2513 


1.1 


1.4 


1.5 


0.9 


1.0 


1.0 


0.6 


0.7 


0.5 


0.2 


0.3 


Alii] 2669 








1.0 


0.8 


0.6 


0.3 


0.3 


0.3 


0.1 


0.2 


He 11 2733 






0.8 


0.8 


0.8 


0.6 


0.6 


0.6 


0.6 


0.3 


0.3 


Fe II 2756 








0.6 




0.4 




0.4 


0.3 




0.1 


Fe II 2774 


0.7 


1.3 


1.0 


1.2 


1.2 


0.7 


0.3 


0.3 








[Mgv] 2783 














0.3 


0.3 


1.0 


1.4 


1.0 


Mgii 2800 


5.3 


5.5 


3.5 


3.1 


8.9 


9.7 


7.9 


4.5 


3.6 


2.1 


1.9 


[Feiv] 2829 




0.6: 


0.7 


0.3 




0.5 


0.6 


0.7 


0.6 


0.4 


0.4 


Oiii 2838 






0.3 


1.8 


2.2 


1.7 


1.9 


1.8 


1.9 


1.0 


0.9 


TTo TT 9SzL7 


A n 


O. i 


1 Q 


n Q 
u.y 


i . 


1 

i .D 


i .u 








n 1 


Fe II 2859 




1.2: 






0.5 


1.1 


0.5 


0.4 


0.3 




0.3 


[Mgv] 2928 


















0.7 


0.5 


0.4 


Mgii 2936 


0.6 


1.0 




0.4 


0.7 


0.6 












[Nev] 2973 














0.5 


0.3 


0.5 


0.4 


0.2 


O III 3025 








0.8 


0.8 


0.8 


0.8 


0.8 


0.9 


0.5 


0.3 


O III 3047 








1.9 


2.2 


1.8 


1.1 


1.8 


1.7 


0.9 


0.8 


Tin 3079 


1.4 


1.6 


1.2 




0.9 


0.4 










0.3 


Oiii 3133 


1.9 


2.1 


6.5 


5.1 


8.3 


9.3 


10.4 


10.0 


9.1 


5.4 


4.7 



There were no measurable emission lines present on JD 2 443 770 and JD 2 449 135 (see also Section 3.1.2). 



2.2 HST Observations 

We carried out iJ^T observations using tiie FOC/96 camera 
on 12 Marcii 1991 (JD 2 448 328). FOC images were taken 
with the F501N ([O III] 5007A) filter and the near-UV ob- 
jective prism (NUVOP) crossed with the F195W fiher. The 
direct image taken with the F501N filter provided measure- 
ments of the star's undispersed position for the analysis of 
the prism observation. 

Unfortunately, the FOC images are not conclusive re- 
garding the point-like nature of RX Pup. Those images with- 
out the prisms in place are, unfortunately, saturated, leaving 
only the very extended halo to hint at the point-like nature 
of the object. The best we can say is that the FOC images 
place an uncertain upper limit to the size of RX Pup of ap- 
proximately 0'.'2. Additional UV and optical images taken 
late in 1993 with the FOC do not provide substantive evi- 
dence of any extended features, down to an upper limit of 
0'.'06. However, the low signal-to-noise in the images and the 
lack of information regarding the actual focus of HST and 
FOC at the time of the observations preclude making any 
definitive conclusions. 

A relatively long exposure time, of 900 s, was required 
for the objective prism image to give sufficiently high signal- 
to-noise across the full spectral range from 1750 - 5500 A. 



The FOC detectors are count-rate limited (unlike CCDs) so 
a 512 X 512 pixel format was used to maintain linearity for 
the star with count rates up to about 0.5 counts/s/pixel. 
This format provides a 11" x 11" field of view, with pixel 
sizes of 0"022 x 0"022. 

The objective prism image taken with PRISM2 -|- 
F195W (NUVOP) provided low resolution spectroscopy of 
RX Pup (Fig. 2). The peak count rates in the objective prism 
image were less than 0.2 counts/s/pixel, insuring photomet- 
ric linearity throughout the entire observed spectral region. 
This observation gave a signal-to-noise ratio ranging from « 
2 at 1800 A to « 14 at 4000 A. 

All of the observations were flat-fielded and geometri- 
cally corrected using the standard FOC pipeline calibrations 
described in detail the FOC Instrument Handbook (Nota et 
al. 1996 henceforth FOCIH). The dispersed images taken 
with the objective prisms are aligned with the visible wave- 
lengths on the lower end and the UV wavelengths at the up- 
per end. The NUVOP image had a resolution of 40A/pixel 
at 5000 A and 0.5 A/pixel at 1700 A extending across the en- 
tire image (512 pixels). 

Seven pixel (0'.'15) wide regions centred on the peak of 
the spectrum were extracted, and the contribution from all 
the pixels summed to provide the raw spectrum. This raw 
spectrum was wavelength calibrated by applying the latest 
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1500 2000 2500 3000 

Wavelength [A] 

Figure 1. lUE low resolution spectra of RX Pup taken in 1978- 
1988. The spectra show the gradual decline of the UV continuum 
and the evolution of emission line fluxes. 




2000 3000 4000 



Wavelength [A] 

Figure 2. HST low resolution spectrum RX Pup taken on 12 
March 1991. 

dispersion relation for each prism as given in the FOCIH. 
The zero point for the conversion was based on the position 
of each star in the F501N image, with the distance along 
the spectrum from that zero-point determining the wave- 
length of the spectrum. Flux calibrating the spectrum re- 
quired some knowledge of the percentage of total flux from 
the source that was contained in the spectrum to recover 
all the flux from the source. This percentage is based on 
the width of the extraction used for the spectrum, which 



was 7 pixels. An empirical relationship for this percentage 
was determined using spectrophotometric standards and the 
method reported in the FOCIH This relationship served as 
the initial estimate for our reductions. With an estimate of 
the total flux contained in the extraction, the raw spectrum 
was then converted from counts to flux units by dividing 
out the FOC and filter sensitivities using software written 
as part of the Space Telescope Science Data Analysis System 
(STSDAS) and described in the FOCIH. 

2.3 Optical spectroscopy 

The optical spectra were taken during 4 observing runs, in 
November 1990 (JD 2 448 200), April 1991 (JD 2 448 352), 
June 1995 (JD 2 449 890) and March 1996 (JD 2 450 167), 
with the 2.15 m telescope of CASLEO at San Juan, Ar- 
gentina. In 1990 and 1991 the BoUer & Chivens Cassegrain 
spectrograph and a Reticon photon-counting system known 
as "Z-Machine" was used. The "Z-Machine" is a double 
aperture detector which permits one to observe simultane- 
ously the sky and the target plus sky. A translation of the 
object from one aperture to the other is made at mid expo- 
sure for each observation, and the combination of these two 
data sets gives a sky-free spectrum. An aperture of 5" x 3" 
was used. The 1200 lines mm~^ grating gives two wavelength 
ranges of about 4400 - 5100A and 5800 - 7200A, with dis- 
persions of 0.5 and 1 A/pixel, respectively. A spectral reso- 
lution of at least 4100 and 2700 in the two spectral regions, 
respectively, is indicated by measurements of the FWHM 
of arc lines. In 1995 and 1996, the spectra were taken with 
the REOSC echelle spectrograph at a resolution of 15 000, 
and recorded with a TEK 1024 x 1024 CCD. The reduction 
of the "Z-Machine" spectra was carried out with the IHAP 
image processing software. Wavelength calibration was per- 
formed using the He-Ar-Ne lamps with reference exposures 
obtained immediately before and after each stellar exposure, 
at the same sky position. Flat-fleld exposures were also done 
and a mean value was used in the reduction procedure. The 
echelle spectra were reduced with standard IRAF packages, 
CCDRED and ONEDSPEC. 

To obtain the flux calibration, standard stars from 
Stone & Baldwin (1983) and Baldwin & Stone (1984) were 
observed each night. A comparison of the spectra of the stan- 
dards suggests that the flux calibration errors are about 15 
per cent for the "Z-Machine" and 20 per cent in the central 
part of each order for the REOSC echelle images, respec- 
tively. 

An additional optical spectrum was obtained in May 
1991 (JD 2 448 385), with the 1.5 m telescope of ESO at 
La Silla, and the B&C spectrograph (used in the range 
4300 - 6900A) equipped with a CCD (512 x 1024 pixels). 
The spectrum was reduced following the IHAP routines, us- 
ing flatfield exposures, HeAr wavelength calibration spectra, 
and standard stars (Hayes & Latham 1975) observed each 
night for flux calibrations. The photometric calibrations are 
accurate to 10-20 per cent. 

Emission line fluxes derived by fitting Gaussian profiles 
are listed in Table 3. These data have uncertainties of ~10- 
15 per cent for strong lines and ~20-30 per cent for weak 
lines or noisy spectra. Examples are shown in Fig. 3. 
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Table 3. Optical emission line fluxes in units of 10 erg cms ^ 



JD 24... 


48200 


48352 


48385 


49890 


50167 


JD 24... 


48200 


48352 


48385 


49890 


50167 


[Fell] 4244 










1.0 


Fell 5317 






1.5 


3.3 


2.7 


C II 4267 










1.3 


[Fell], Fell 5362 






0.7 


2.2 


1.7 


[Fell] 4287 










1.8 


[Fell], Fell 5413 






1.1 


1.2 


1.5 


H7 4340 






5.5 




6.5 


Fe II 5535 






0.6 


2.3 


2.2 


[Oiii] 4363 






2.3 




0.7 


[Oi] 5577 






0.5 


0.5: 


1.0: 


He I, Fe 11 4387 






0.5 




0.6 


[Nil] 5755 






8.0 


1.0 


0.7 


[Fell] 4414 


1.6 


0.8 


1.4 




1.0 


He I 5876 


61.0 


4.9 


4.2 


0.1: 


0.4 


Fell 4417 




1.0 






0.7 


Nai 5890 




3.7 


3.1 


abs 


abs 


[Fell] 4452 




0.4 


0.7 




1.1 


Nai 5896 




2.6 


2.8 


abs 


abs 


[Fell] 4458 




0.5 


0.6 




0.4 


Fe II 5991 






0.6 


3.0 


3.9 


He I, Fe II 4472 


3.1 


1.1 


1.1 




0.3 


Fe II 6084 






0.4 


2.3 


1.9 


[Fell], Fell 4491 




0.6 


0.8 




0.6 


[Fevii], [Cav] 6086 


86.0 










Hen 4541 


1.1 










Fell 6149 






0.8 


2.9 


3.4 


Fe II 4584 






1.3 




2.8 


Fe II 6240 






0.7 


5.5 


3.8 


[Feiii] 4607 


0.6 










Fe II 6248 






0.8 


3.5 


2.9 


N III 4634 


1.7 










[Oi] 6300 


31.8 


29.0 


25.1 


16.9 


13.0 


N III 4641 


4.0 




1.1 






[Siii] 6310 


29.0 










[Feiii] 4658 


3.0 


1.1 


0.8 






[Oi] 6364 


11.8 


10.1 


8.2 


3.9 


4.1 


He II 4686 


39.1 


1.1 


1.1 






Fe II 6433 






1.0 


7.1 


6.3 


[Feiii] 4702 


2.2 




0.2 






[Arv] 6435 


7.5 










[Ne iv] , He I 4714 


3.8 




0.5 






Fe II 6456 






1.2 


4.4 


4.8 


[Neiv] 4725 


4.7 




0.9 






Fell 6516 






1.4 


5.3 


7.1 


[Fell] 4815 


0.6 


0.7 


0.7 




1.1 


[Nil] 6548 




28.1 


20.6 


30.0 




H/3 4861 


110.1 


19.2 


18.8 


29.0 


20.0 


He 6563 


1790.0 


173.6 


145.8 


700.0 


630.0 


[Fell] 4890 


0.8 


0.5 


1.0 


0.8 


0.6 


[Nil] 6584 


79.1 


78.2 


64.4 


60.0 


53.4 


[Fevii] 4894 


1.1 










He I 6678 


18.6 


2.2 


2.0 




1.0 


[Fell] 4905 


1.7 






0.5 


0.4 


[S 11] 6717 


3.8 


2.8 


2.2 


1.1 




He I 4922 


2.0 










[Sii] 6731 


5.8 


5.7 


4.5 


2.7 




Fe II 4924 




1.5 


1.6 


4.0 


5.5 


[Foiv], [Kiv] 6793 


2.1 










[Fevii] 4942 
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2.4 Photometry 

JHKL broad-band photometry (1.25, 1.65, 2.2, 3.45 ^m) 
was obtained with the Mk II infrared photometer on the 
0.75 m telescope at SAAO, Sutherland (see Carter 1990 for 
details about the system). Some of the early data have been 
published (W83; Whitelock et al. 1984, hereafter W84), al- 
though the values given here (Table 4) might be slightly 
different from those published as they have been revised fol- 
lowing the work by Carter (1990). The uncertainty on indi- 
vidual measurements is less than 0.03 mag in JHK, and less 
than 0.05 mag at L. Less precise observations are marked 
with a colon. 

When available, the FES counts were converted to V 
magnitudes using the time dependent correction from Fire- 
man & Imhoff (1989) and the colour dependent equations 



from Imhoff & Wasatonic (1986). The resulting magnitudes 
given in Table 1 have accuracies of 0.1 mag. We have also col- 
lected all published optical photometric data (Eggen 1973; 
Feast et al. 1977; Klutz 1979; W83, W84; Munari et al. 1992; 
IS94) as well as visual magnitude estimates reported on reg- 
ular basis by the Variable Star Section of the Royal Astro- 
nomical Society of New Zealand (RASNZ). Fig. 4 illustrates 
the V/nivis light curve together with the IR data. 



3 ANALYSIS 
3.1 Variability 
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Table 4. Infrared photometry of RX Pup 
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Figure 3. Optical spectra of RX Pup. 



3.1.1 A brief survey of observations prior to 1969 

Although the first observational records of RX Pup date 
back to the end of XIX century, the star has been generally 
neglected by most observers, and only sketchy information 
on its spectroscopic and photometric variations prior to the 
early 1970s is available. 

RX Pup was classed among "stars having peculiar spec- 
tra" by Fleming (1912). The spectrum taken on November 
24, 1894 resembled that of 77 Car, showing bright H/?, H7, RS 
and H(" (Pickering 1897). Another spectrum taken on Jan- 
uary 8, 1897, showed in addition to the strong H Balmer 
lines, bright Hell A4686, [Oiii] A4363 and [Neiii] A3869 
(Cannon 1916). Cannon also noted that the spectral lines 
(although not their widths) resembled those of novae in the 
nebular phase. Extensive photometry obtained in 1890-1924 
(Pickering 1914; Yamamoto 1924) showed the star bright- 
ened from TTipg = 12 in 1890 to mpg = 11.13 in 1894, then 
declined to a very deep minimum (m-pg = 14.1) in 1904-5. 
After recovery c.1909, the photographic magnitude showed 
only small fluctuations between 11.6 and 12.3 mag. 

The next spectra of RX Pup, obtained in January 1941, 
revealed a continuum confined to the yellow and red region, 
and very strong high excitation emission lines of [Fevil], 
[Ne v] , [Fe vi] and [Ca vil] , in addition to strong H I Balmer 
series. Hell, [Nelll], and [Om] emission lines (Swings & 
Struve 1941). Swings & Struve also pointed out the spectro- 
scopic similarity of RX Pup to the symbiotic stars CI Cyg 
and AX Per, although they did not find definite evidence for 
a late- type component in RX Pup. 

The spectra taken between 1949 and 1951 showed very 
faint (or absent) continuum, moderately sharp Ha, as well 
as [0 1] AA6300, 6364 and other blended emissions (most of 
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Figure 4. Optical and IR light curves of RX Pup in 1968-1998. 
Dots: visual observations from RASNZ, open circles: published V 
mags, crosses: FES mags, square with a down arrow: February 
1968 visual mag limit. 



them probably Fe 11 lines) between A 6300 A and Hq (Henize 
1976). A similar low ionisation spectrum, with no Hell 
A4686, [O III] A5007 very weak relative to H/3 and a very faint 
continuum, was observed on February 26, 1968 by Sanduleak 
& Stevenson (1973) who also estimated the visual and pho- 
tographic magnitudes at rrivis ^ 12.5 and rripg = 13.5 ± 0.5, 
respectively. 



3.1.2 The spectral transformations in 1969-1996 

Between 1968 and 1969 the optical brightness of RX Pup 
increased by a few magnitudes (Fig. 4). Between Decem- 
ber 1969 and February 1970, the star reached V ~ 9.5 — 9.3 
(Eggen 1973). Its colours B-V ~ 1.3-1.2, U-B ~ 0.5-0.4, 
and R — I ~ 0.8, corrected for a reasonable amount of the in- 
terstellar reddening, Eb-v ~ 0.6 (Section 3.2), indicated the 
presence of a bright blue continuum. Optical spectroscopy 
obtained in February 1972 and March 1975 (Swings & Klutz 
1976; also Webster & Allen 1975) revealed a bright late A- or 
early F-type continuum with strong H I as well as O I, [O l]. 
Can, [Can], and numerous Fell and [Fen] emission lines. 
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Figure 5. Variation of UV continuum magnitudes in 1978-1993. 



RX Pup maintained this low excitation emission spectrum 
together with a strong blue continuum until March 1979 
(Klutz, Simonetto & Swings 1978; Klutz & Swings 1981). 
The optical brightness reached a maximum in 1977-78, with 
V ~ 8.5, B - V 0.9, and [/ - B ~ 0.1 - 0.5, the op- 
tical spectrum resembled that of a Be or shell star (Feast 
et al. 1977; Klutz 1979). The first lUE spectrum taken in 
September 1978 (Fig. 1) shows a UV continuum like that 
of a Be/shell star. Similar optical/UV spectra have been 
observed in CH Cyg and other symbiotic stars during the 
visual maximum of their outbursts (Mikolajewska, Selvelli 
& Hack 1988; Kenyon 1986). 

The H I Balmer lines presented double P Cygni profiles 
(V << R) in 1972, which evolved into a complex P Cygni 
structure in 1975, with one emission and a series of absorp- 
tion components ranging from -1100 to -230 kms~^ (Swings 
& Klutz 1976). Such complex structure was also observed in 
1976, although higher members (such as H7 and HJ) showed 
very little absorption. In 1978 the Balmer lines became fairly 
stable emission lines with a sharp P Cygni profiles and a few 
fainter absorptions with positive velocities (Klutz 1979). 

In December 1979 RX Pup returned to the high excita- 
tion conditions exhibited 40 years earlier; the appearance of 
strong and broad emission lines of Hell, Nlll and [Qui] in 
the optical range was accompanied by decline in the UBV 
fiuxes (Klutz & Swings 1981; see Fig. 4). The Balmer contin- 
uum strengthened, and the intensities of the A4640 feature 
and He 11 A4686 emission resembled those observed in WN7- 
8 stars. The 1979 spectrum, however, showed neither [Fevi] 
nor [Fevil]. The next optical spectrum taken in April 1983 
revealed faint [Fcvi] and [Fevil] emission lines (W83), with 
a flux ratio [Fe vil] A6086/He I A5876 ~ 0.2. The intensity 
of the [Fevil] was growing, and in the 1984/89/90 spectra 
(AW88; de Freitas Pacheco & Costa 1992, hereafter FPC92; 
Table 3) the [Fevii] A6086/He I A5876 flux ratio increased to 
0.3/1.7/1.4, respectively. The changes in the optical spec- 
trum observed between 1979 and 1989 clearly indicate a 
graduaUy increasing degree of ionisation, similar to that re- 
ported for symbiotic novae (e.g. Kenyon 1986). 

The evolution of the blue continuum and emission line 
spectrum in 1978-1993 can be traced best using lUE data. 
Fig. 1 illustrates the changes in the UV. To establish the 
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Figure 6. Evolution of optical and UV emission line fluxes in 
1979-1996. The Ha and [O III] A5007 fluxes published by W83, van 
Winckel, Duerbeck & Schwarz (1993) and IS94 are also plotted. 



continuum behaviour, we measured average fluxes, F\, in 
20 A bins free from emission lines and derived magnitudes 
using: m\ = — 2.51og_FA — 21.1. The continuum declines 
steadily from 1978 to 1991, in a similar way to the visual 
light (Figs 5 & 4); the UV continuum seems to fade less 
at short wavelengths (AA 1500A) than at longer wave- 
lengths (AA 25OOA). The continuum decline is matched 
by an increase in flux from most of the emission lines (Fig. 
6). In particular, the Hell A1640 flux increased by a fac- 
tor of ~ 5 between June 1979 and AprU 1980. The Hell 
flux continued a slower increase until June 1981 and from 
then to the observation of May 1986 it remained at roughly 
constant level. Similar behaviour was seen amongst other 
medium excitation lines, such as Civ A1550, Niv] A1487, 
Si IV A1394 and Oiv] A1403. Nv A1240 appeared for first 
time in June 1981, while [Mgv] A2783 appeared in Octo- 
ber 1983 (Fig. 7); these lines showed maximum fiuxes in the 
spectra taken during 1986 and 1988, respectively. Generally, 
the trends observed in the UV continuum and emission line 
fiuxes show an evolution towards higher stages of ionisation, 
similar to that observed in symbiotic novae (Nussbaumer & 
Vogel 1990; Miirset & Nussbaumer 1994). 

From an analysis of high resolution lUE spectra taken 
between September 1980 and March 1984, K85 (also K82) 
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Figure 7. Appearance and evolution of [Mgv]. The high resolu- 
tion spectra were taken on (from bottom to top): 22 March 1982, 
30 October 1983, 9 May 1986 and 25 October 1988, respectively, 
show the appearance and increase of [Mgv] A2783. 



found that many of the strong emission lines exhibit mul- 
tiple component structure. Our analysis, however, does not 
confirm some of their identifications. For example, the red- 
shifted component of Hell A1640 as well as the blue-shifted 
component of Nlll] A1749.7 (Fig. lb and Fig. Id of K85) 
are misidentified Ol] A1641.4 and Nlll] A1748.6, respec- 
tively. The intensity of the A1641 feature is strongly cor- 
related with the O I A 1304 line flux measured on low reso- 
lution spectra; this feature is absent on all high resolution 
spectra taken in 1983-1988 when the Ol A1304 was very 
faint or absent. Moreover, profiles of all semiforbidden emis- 
sion lines detected on the high resolution spectra show more 
or less the same velocity structure, with the width at the 
foot of ~ 250kms~^ centred roughly at the stellar velocity 
(~ lOkms"^; Wallerstein 1986). The only exception would 
be N III] A1749.7 if the K85 interpretation were accepted. We 
also disagree with K85 about the close similarity between the 
Civ and Hen profiles. The Civ AA1548.2, 1550.8 doublet 
emission profile is truncated by a narrow P Cygni structure 
at ~ — 90kms~^, and it exhibits rather complex structure, 
while He ll has a considerably simpler profile which is almost 
symmetric around the stellar velocity. 

The C IV doublet shows unusual intensities, the 
7(A1548.2)//(A1550.8) ratio had an observed value of ~ 0.6 



in 1980-82, then in 1983 its value increased to ~ 0.8. In 
1987 the Civ doublet ratio reached the optical thick limit 
of unity (see also Michalitsianos et al. 1988), while in 1988 
/(A1548.2)//(A1550.8) ~ 1.1. The Civ ratio seems to in- 
crease following the general enhancement of high-excitation 
emission line fiuxes (Table 2; Fig. 6). Similar behaviour 
was reported in the symbiotic nova HM Sge (Miiller & 
Nussbaumer 1985). Michalitsianos, Kafatos & Meier (1992) 
demonstrated that the combined fiuxes of Fell Bowen- 
pumped lines can account for an appreciable portion of the 
flux deficit in the C iv A1548.2 line in 1980-83, when the C iv 
doublet ratio is less than the optically thick limit. Their in- 
terpretation is consistent with the presence of O l] A1641.3 in 
1980-83 indicating the emitting region was optically thick. 

By early 1991, the high excitation emission lines had 
practically disappeared. The HST spectrum taken in March 
1991 (Fig. 3) revealed a very faint UV continuum without 
strong emission lines. In particular, the intercombination 
Silll]A1892 and Clll]A1909 stiU strong on the SWP spec- 
trum of May 1989 had disappeared, as had the O iii A3133 
emission which was very strong in the 1980s. The only emis- 
sion detected in the UV range was a weak Mg ll A2800. The 
optical continuum was also faint, consistent with the deep 
minimum shown by the V/nivis fight curve at this time (Fig. 
4). The lack of measurable short- wavelength continuum or 
emission lines (except Mg ll A2800 and very faint Fe ll A2756) 
was confirmed by the lUE spectra taken a few days after the 
HST data. The optical spectra taken in April and May 1991, 
also showed a significant decline in the continuum and emis- 
sion line fluxes (Fig. 3 and Fig. 6). The only emission lines 
that remained practically unchanged were low excitation for- 
bidden lines [Ol], [Nil] and [Fell]. The 1991 spectrum was 
also characterised by the appearance of NalD in emission 
(Fig. 8). The decrease in the Hell/H/3 ratio from ~ 0.4 in 
the 1990s to ~ 0.06 in 1991 suggests a general decline of the 
source of ionisation. 

Between 1992 and 1993, the optical magnitudes recov- 
ered to about the same level as in the late 1980s (Fig. 4), 
the degree of excitation, however, continued to decrease. 
He II A4686 was absent on optical spectra taken in 1995 and 
1996, and the fiuxes in practically all of the emission lines 
were significantly reduced. The only lines that remained rel- 
atively strong were the Balmer H I lines (Ha brightened by 
a factor of ~ 4 since 1991; Fig. 6), although weak low- 
excitation lines from Fell and forbidden transitions were 
also visible. The Balmer decrement was unusually steep, 
Ha:H/9 ~ 30, and both Ha and H/9 had deep central ab- 
sorptions (Fig. 9). 

RX Pup also exhibits interesting variability at radio 
wavelengths. In the mid-1970s, when detected for the first 
time, the radio emission was characterised by a flat, opti- 
cally thin spectrum. By 1985 it had changed to a steeply 
rising spectrum (Seaquist & Taylor 1987). In the 1990s, it 
flattened at all wavelengths with the previous level of flux 
density maintained only at the longest wavelength (IS94). 
The transition from a steep to a fiat radio spectrum seems to 
have occurred between February and December 1991 (IS94). 
The 3.5 cm flux dropped from Sg.sGHz ~ 69.7 ± 3.5 mjy 
in February 1991 (Seaquist, Krogulec & Taylor 1993) to 
Sg.gGHz = 11.3 ± 0.3 mJy in July 1992 (ATCA; Lim 1993) 
and was apparently related to changes in the UV and optical 
spectra. 
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Figure 8. Evolution of Nal Di,D2 lines. 

In contrast to the optical brightness, the infrared 
JHKL light curves are dominated by periodic variations 
consistent with those expected for a Mira variable, and 
a long term, wavelength dependent, decrease of the mean 
brightness (Fig. 4). A Fourier analysis of the data in Table 
4 confirms the pulsation period of 578 day found by Feast 
et al. (1977). The following ephemeris gives the phases of 
minima in the J band: 

Min( J) = 2 440 214 + 578 X (1) 

The long term trends were attributed to dust obscuration 
possibly related to the orbital motion (W84; Whitelock 
1987). The IR fading is strongest in the J band, and de- 
creases in amplitude with increasing wavelength. We es- 
timate the mean fading rates in 1979-1991, AJ ~ 0.2, 
AJ/" ~ 0.12, and AiC ~ 0.08magyr"\ respectively. The 
fading started in late 1979, following the decline in vi- 
sual brightness and appearance of the high-excitation emis- 
sion line spectrum. The faintest infrared flux was measured 
around 1991, after which the JHKL magnitudes recovered 
to the values observed in 1982-1983. 

Several authors claimed that the visibility of molecu- 
lar bands from the Mira is limited to the periods of low 
excitation (IS94, and references therein). In fact, although 
RX Pup was undoubtedly in the low excitation phase in the 
1970s, red and blue spectra taken in 1972 show no trace of 
absorption bands (Swings & Klutz 1976). Similarly, Swings 
& Klutz failed to detect any TiO bands at A 86OOA on 
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Figure 9. Hi Balmer lines in 1995 and 1996 (thick and thin 
line, respectively) exhibited strong central absorptions which were 
probably responsible for the very steep decrement. The Ha profile 
observed in 1988 (dotted line) by van Winckel et al. (1993) is 
shown for comparison. 

their spectra taken in March 1975. Strong TiO bands in the 
AA7589 - 8432 range as well as VO A7900 were detected, 
for the first time, on November 30, 1980 (Andrillat 1982), 
after RX Pup entered a high excitation phase (!). Although 
Schulte-Ladbeck (1988) did not find any TiO at A < 8OOOA, 
the TiO as well as VO bands were visible beyond A ~ 8OOOA 
in April 1984 (AW88). Similarly, there were no detectable 
absorption bands on our red spectra taken in 1990 and 1991 
(A ^ 72OOA), while the absorption bands beyond A ~ 7200A 
were present on spectra taken in 1993 by IS94, and on our 
spectra taken in 1995 and 1996. In all cases, the absorp- 
tion bands of the Mira were detected beyond A ~ 7000A. In 
particular, the TiO bands at A6I8OA and A7100A, strong in 
most symbiotic systems (Kenyon & Fernandez-Castro 1987), 
were never detected in RX Pup, including our well exposed 
spectra of 1995 and 1996. The only positive detections were 
made in the long wavelength range, A ^ 72OOA, when the 
visual magnitude of RX Pup was V ~ rrivis ^ 10.5 and near 
maxima of the Mira's pulsations. The permanent absence 
of the molecular absorption bands in the optical range is 
consistent with the fact that the Mira pulsations were never 
visible in that range (Fig. 4). Similar behaviour is shown by 
other symbiotic Miras. In particular, Mira pulsations were 
never observed in the optical light of HM Sge, even during 
quiescence (Yudin et al. 1994). 
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Figure 10. Spectral development of the 1969-1990 outburst of 
RX Pup. Solid, dashed and dot-dashed lines indicate emission and 
absorption features, and complex P Cyg profiles, respectively. 



The evolution of the B — V and U — B colour indexes 
in RX Pup also suggest that the contribution of the Mira 
to the optical light was never significant. The B ~V index 
had roughly constant value ~ 0.9 in 1975-1990, which in- 
creased to ~ 1.7 in 1993, while the U — B index increased 
from ~ 0.1 in 1975 to ~ 0.5 in 1979, then decreased to 
~ —0.5 in 1989-1990, and increased again to ~ 0.4 — 0.5 
in 1993 (Feast et al. 1977; Klutz 1979; Munari et al. 1992; 
IS94). AW88 found the 1984 optical continuum resembled 
reddened gaseous emission. We estimate contribution of the 
emission lines to the observed broad band B and V magni- 
tudes, ~ 1/0.25/0.05 and AS - 1.8/0.45/0.05 mag in 
1990/91/95, respectively. During the 1993 observation the 
Mira component was near maximum of its pulsation vari- 
ability (0 ~ 0.9). Given the reddening towards RX Pup, 
Eb-v 0.6 (Section 3.2), and that near maximum light 
long period Miras have B — V ^ 1.6 and U — B 1 (Smak 
1964), the observed values of (7 - B and B - V for RX Pup 
are similar to those found in other symbiotic systems, and 
indicate the hot component dominates at these wavelengths. 

Fig. 10 summarises the spectral development of RX Pup 
during the last three decades. The parallel between RX Pup 
and the symbiotic novae AG Peg (Gallagher et al. 1979; 
Kenyon et al. 1993), V1329 Cyg (Kenyon 1986, and refer- 
ences therein), HM Sge (Nussbaumer & Vogel 1990) as well 
as the symbiotic recurrent nova RS Oph (Kenyon 1986) is 
apparent, although these other systems seem to be evolving 
on very different time scales. 



3.2 Reddening and distance 

Before deriving the physical characteristics of the RX Pup 
binary, we need a good estimate for the extinction and dis- 
tance. We base our reddening estimates on the UV contin- 
uum, UV and optical emission line ratios, IR photometry, 
and measurements of the interstellar features. We also as- 
sume the reddening to be represented by a standard inter- 
stellar extinction curve (e.g. Seaton 1979). 



3.2.1 Interstellar absorption features 

RX Pup lies in field 99(256/-3.0) of Neckel & Klare (1980), 
for which they estimate the interstellar reddening to be 
Eb-v ~ 0.6 at d ~ 1.3 — 2 kpc. For such an E_g-y we 
would expect measurable interstellar features. 

In fact. Klutz, Simonetto & Swings (1978) reported in- 
terstellar Nal Di,D2 absorption lines as well as broad ab- 
sorption bands at A ~ 5780 and 5797 A. The total equiv- 
alent width, EW(Di + D2) ~ O.43A, measured by Klutz 
et al. implies ¥^b-v ^ 0.1 (Sembach, Danks & Savage 
1993; Munari & Zwitter 1997), while EW(A5780) = 0.45A 
and EW(A 5797) = 0.52A combined with calibrations of 
Josafatsson & Snow (1987) and Herbig (1993) indicate a 
much higher E_B-y ~ 0.8 — 1.0. The latter value is com- 
parable to our estimates based on emission line ratios (see 
below), but higher than the interstellar reddening predicted 
by Neckel & Klare (1980). 

The absorption lines of Nal Di,D2 are also clearly vis- 
ible in our optical spectra taken in 1995 and 1996 (Fig. 8). 
The profiles are, however, very complex, with faint emission 
features truncated by very strong absorption components, 
EW(Di) - EW(D2) - lA, at ~ -2 km/s, and weaker 
components, EW(Di) ~ 0.3 - 0.4 and EW(D2) ~ O.lA, at 
~ —180 km/s, on both spectra. Traces of both absorption 
components are also visible in the emission line profiles of 
Nal on the spectrum taken in April 1991. The ESO spec- 
trum (May 1991) has insufficient resolution to show them. In 
addition, the 1995 spectrum shows strong Kl AA7665, 7699 
absorption lines [EW(7665) - 1.3A, EW(7699) ~ 0.8A] at 
~ 1 km/s, very close to the radial velocity of the stronger 
components of the Nal D lines. It seems unlikely that these 
strong absorption features are of interstellar origin. The ra- 
dial velocities of the K I absorption lines and of the stronger 
components of the Nal lines are ~ —16 km/s and —19 km/s, 
respectively, with respect to the local standard of rest (LSR) 
and they are blueshifted with respect to the velocity pre- 
dicted by galactic rotation law (1; ~ 5 — 10 km/s for a rea- 
sonable range d ~ 1 — 2 kpc). The Kl line is somewhat 
asymmetric, while the Nal absorption lines seems to consist 
of 2 components. The strong absorption components of the 
Nal doublet as well as the Kl absorption lines are probably 
superpositions of the interstellar and circumstellar/stellar 
components. 

Adopting the calibration of Bromage & Nandy (1973), 
we estimate from the combined equivalent widths of both 
Nal D lines (EW(Di -I- D2)that ~ 2A) Es-v ^ 0.6. Other 
calibrations give similar results, for example Munari & Zwit- 
ter (1997) find EW(Di) - lA for 'Eb-v ^ 0.5 and stars with 
multi-component interstellar lines. However, the K I absorp- 
tion lines are too strong for any plausible E_b-v (e.g. Chaffee 
& White 1982; Munari & Zwitter 1997). 

The spectra taken in 1995 and 1996 also show the dif- 
fuse interstellar bands at A ~ 5780, 5797 and 6284 A. Their 
equivalent widths, EW(A5780) ^ 0.2lA, EW(A5797) « 
O.I2A, and EW(A6284) 7» 0.22A, are consistent with 
¥.B~v ~ 0.5 (Herbig 1993; Josafatsson & Snow 1987; Snow, 
York & Welty 1977). 
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3.2.2 UV continuum 

An independent reddening estimate towards the hot com- 
ponent of RX Pup can be derived from the intensity of the 
A2200 band. Unfortunately, practically all lUE spectra have 
very poorly exposed continua. The only spectra that can 
be used are SWP 2684 and LWR 2395, taken near the vi- 
sual maximum in 1978. These show characteristic features 
of the iron curtain absorption and a strong absorption band 
around A2200A. We compared the combined SWP -f LWR 
spectrum with the lUE spectra of CH Cyg taken during its 
1981-1984 visual maximum (Mikolajewska et al. 1988), and 
estimated 'Eb-v ^ 0.6. 

The very strong absorption feature near A2200A is also 
visible on the HST spectrum taken in May 1991. Correcting 
this spectrum for Eb-v ~ 0.8 ± 0.1 entirely removes the 
feature, and the dereddened spectrum resembles that of CH 
Cyg observed during the 1985 minimum (Mikolajewska et 
al. 1988). Given the uncertainty in the shape of the hot 
component UV continuum, these estimates generally agree 
with those derived from the emission lines. 



3.2.3 Emission lines 

Both the [Nev] AA 1575, 2975 lines, originate from the 
same upper level, as do the [Fevil] AA 5721, 6086 lines. 
Thus their intensity ratio should depend only on the tran- 
sition probabilities, [Ne v] /(A1575)//(A2975) = 1.5, and 
[Fevil]/(A5721)/7(A6086) = 0.65, respectively. The average 
[Nev] ratio, ^ 0.6 measured on the lUE spectra taken in 
1983-1988, implies Fjb-v = 0.4 ± 0.4. The high uncertainty 
being due to the faintness of the [Nev] lines. The [Fe vil] line 
ratio observed in 1983 by W84 indicates Eb-v = 1.4 ± 0.5, 
and data from FPC92 give Es_v = 1.9 ± 0.5 in 1989. 
These values should be however considered as upper lim- 
its for Es_v because the [Fevil] A6086 line can be strongly 
blended with the [Cav] A6086 line. 

A better reddening estimate can be made from the 
[Neiv] A1602/A4725 ratio (e.g. Schmid & Schild 1990a). The 
intensity of the ultraviolet [Neiv] A 1602 emission line was 
practically constant in 1980-89. We compared its average 
flux with that in the optical [Ne iv] A4725 line ( W84) , and 
estimate Eb-v = 0.9 ± 0.1. 

The [On] /(A7325)/7(A2471) ratio has a practically 
constant value of 1.27 for electron temperatures. To, be- 
tween 5 000 and 20 000 K and densities, rie, from 10 to 
10^ cm"^ (Canto et al. 1980). The A 2471 line is absent 
on the lUE spectra taken before 1983, while a very weak 
emission appears on later spectra. The optical [O ll] A7325 
doublet was measured on our 1995 and 1996 spectra, and 
it seems to be also present on the 1993 spectrum (IS94). 
The A 7325 emission was also present on the 1984 spectrum 
published by AW88, and on the 1989 spectrum measured 
by FPC92. Our estimate of /(A7325)//(A2471) ^ 40 in 
1984 and 1989, is consistent with Eb-v ^ 0.7. Similarly, 
/(A7325)//(A2471) ^ 13 in 1993-1996, suggests Eb-v ^ 
0.5. 

Extinction parameters can also be estimated from opti- 
cal H I, He I and He ll recombination lines. The ratios of He ll 
A1640, A2733 and A3203 lines observed in 1982-1988 are con- 
sistent with case B recombination and Eb-v = 0.7 ± 0.1. 
The Hell A1640 emission line flux was practically constant 



in 1982-1984. Combining the average flux in the UV line 
with that in He ll A4686 measured in 1983 by W84, we esti- 
mate I(A1640)//(A4686) « 0.3 on JD 2 445 459, which indi- 
cates Eb-v ~ 0.8. Finally, the published Hell A4686/A5411 
ratios indicate Eb-v ~ 1-4 (W84), and ~ 1.2 (FPC92), re- 
spectively. These values should be, however, considered as 
upper limits because the A5411 line is blended with Feii 
A5413, and [Feiv] A5426. 

The H I Balmer lines are reliable reddening indicators 
only for negligible self-absorption in the lower series mem- 
bers (case B recombination, e.g. Cox & Mathews 1969; Net- 
zer 1975). For most D-type symbiotic systems reddening val- 
ues estimated assuming case B recombination agree with 
those derived using other methods (Mikolajewska, Acker & 
Stenholm 1997). In the case of practically all S-type systems, 
and some D-types, there are significant departures from case 
B, and the reddening- free Hci/H/3 ratio is ~ 5 — 10 (e.g. 
Mikolajewska & Kenyon 1992a; Mikolajewska et al. 1995; 
Proga et al. 1996; Mikolajewska et al. 1997) due to self- 
absorption. In particular, Ha/H/3 ^ 5 — 10 is found at the 
high electron densities, no ~ 10^ — lO^^cm"^ (Drake & Ulrich 
1980) which characterise the Balmer line formation region 
(see Section 3.4). 

The values of Ha/H/3, H7/H/3 and H(5/H/3 from W84 
and FPC92 are inconsistent with case B recombination for 
any reddening. If the Ha optical depth is thq ~ 5 — 10, 
the lines ratios imply Eb-v ~ 1- The Hi ratios measured 
in 1991 (JD 2 448 385; Table 3) are consistent with case B 
recombination, and Eb-v ~ 1-0 ± 0.1. Finally, our observa- 
tions in 1995 and 1996 show a very steep Balmer decrement, 
Ha/H/3 — 24 and 31, respectively. Similarly, IS94 measured 
Hq/H/3 = 31 in June 1993. IS94 also derived Eb-v ~ 0.4, 
comparing the flux in H/3 with the radio fiux at 5 GHz ob- 
served at the same epoch. However, their estimate assumes 
case B recombination, and is entirely inconsistent with the 
Eb-v ~ 2.3 suggested by the Ha/H/9 ratio. The Hi emis- 
sion line profiles observed in 1995 and 1996 are complex with 
central reversals (see also Fig. 9). The Ha, H/3 and H7 lines 
are apparently affected by self-absorption. The higher mem- 
bers are generally less affected. In particular, the H7/H/3 
ratio differs from the case B value by less than 10 per cent 
as long as thq ^ 100 (Netzer 1975; Proga et al. 1996). As- 
suming that the bulk of H I Balmer emission originates from 
the same region as the optical He I and UV intercombination 
lines, we estimate rua ^ 20 in the 1980s and 1995-96, and 
^ 0.5 in 1991. The observed H7/H/3 ratio was practically 
constant in 1983-1996 within the observation accuracy of 5- 
10 per cent, and its mean value 0.32±0.01 is consistent with 
Eb-v = 0.8 ± 0.3. 

The calculations of Proga, Mikolajewska & Kenyon 
(1994) demonstrated that the Hei 7(A6678)/7(A5876) emis- 
sion line ratio distinguishes between S- and D-type systems 
in that it is ~ 0.25 for D-types and ^ 0.5 for S-types. The 
average He I 7(A6678)//(A5876) = 0.38 ± 0.03 for RX Pup 
indicates Eb-v = 1.0 ± 0.3. 

The emission line ratios thus generally indicate the red- 
dening of Eb-v ~ 0.8 to 1.0 during the period 1980-1991. 

3.2.4 Mira component 

It is clear (see Section 3.1) that the Mira component made 
a negligible contribution to the optical light, at least during 
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the last three decades. Therefore, estimates of the redden- 
ing toward the Mira must be based on an analysis of the 
infrared. Since the IR light curves (Fig. 4) shows large vari- 
ations in the average light and colour of the Mira, we made 
independent estimates for 4 different epochs corresponding 
to different stages of the hot component (see also Section 
3.3): 

epoch 1: JD 2 442 767-4 000: the optical maximum and low 
excitation phase. 

epoch II: JD 2 445 000-6 600: constant hot component lu- 
minosity and high excitation phase. 

epoch III: JD 2 447 500-8 500: optical and IR minimum. 

epoch IV: JD J> 2 448 500: quiescence and low ionisation. 

The colours of Miras have been discussed by Feast et al. 
(1989) for the LMC, Whitelock et al. (1994) for the South 
Galactic Cap (SGC) and Glass et al. (1995) for the BulgeQ 
These results suggest that the colours of the Miras depend to 
some extent on the environment in which they occur, prob- 
ably as a result of the different metallicities of the various 
environments. However, (J — K) appears to be less sensitive 
to metallicity variations than are other colours. It also pro- 
vides a longer baseline for reddening estimates than {J ~ H) 
or {H — K) and avoids any contribution from dust emis- 
sion which may effect the L observations. We therefore use 
(J — K) in this analysis in preference to other IR colours. 

The reddening towards the Mira was estimated using 
the (J — K) period-colour relation derived by Whitelock et 
al. (1998 in preparation) for Miras in the solar neighbour- 
hood: 

(J-_ft:)o = -0.78-f 0.871ogP (2) 

Thus for RX Pup, (J — A")o = 1.63. The values derived are 
listed in table 5. 

W83 suggested that there may have been a hot source 
contributing to the IR during what we are now calling phase 
I. We also note that the Mira pulsation amplitude was low 
(AJ - 1, AJs: ~ 0.8 in epoch I, and AJ ~ 1.4, AK ~ 1 
during remaining epochs) at this stage - consistent with the 
additional contribution from a warm non-pulsating source. 
We therefore deduce that during phase I, the hot component 
contributed significantly to the J flux. The value of E_g-v 
tabulated for this phase is therefore likely to be an under- 
estimate. A similar effect may have been observed in HM 
Sge, where the amplitude of Mira pulsations decreased from 
A J ~ 1.2 at the visual maximum during the hot component 
outburst to A J ~ 2 mag during the late outburst decline. 

Unfortunately, almost all well-studied long period Mi- 
ras (P ^ 400"*) experience circumstellar reddening, although 
some have much more than others (e.g. Whitelock et al. 
1994; Glass et al. 1995). Three Miras in the SGC have pe- 
riods similar to RX Pup, they are AW Psc, AFGL 5592 and 
NSV 14540 (Whitelock et al. 1994). The colours of these 
stars ((J — K) = 1.9 to 2.4) may be presumed to be un- 
affected by interstellar reddening. Thus, if we assume in- 
terstellar (or interstellar plus circumbinary) reddening pro- 
duces Ej_K ^ 0.4 (corresponding to "Eb-v = 0.8), it ap- 

§ Note that these papers present data on two different photo- 
metric systems, thus transformations must be made before stellar 
colours can be compared 



Table 5. Reddening towards the Mira component of RX Pup 



Epoch 


K 


J -K 
(mag) 




I 


2.37 


2.13 


> 1.0 


II 


2.90 


2.72 


2.3 


III 


3.68 


3.24 


3.3 


IV 


3.12 


2.58 


2.0 



pears that the colours of RX Pup were quite normal during 
phases II and IV. During phase III, however, the Mira was 
considerably redder. 

We also find that the IR continuum energy distribution 
of RX Pup observed by IRAS (during epoch II) is practically 
identical to the IR continua of the long period SGC Miras 
(AW Psc, AFGL 5592 and NSV 14540) mentioned above. 
This supports the interpretation that, during phases II and 
IV, the Mira in RX Pup was a normal high-mass-loss AGB 
star. Whitelock et al. (1994) also show correlations between 
mass-loss rates and K — [12] colour, pulsation period, and 
the amplitude of the if-light curve for their SGC sample. 
The observed K — [12] « 4.6 (the K magnitude has been 
corrected for E_b_v ~ 1) for RX Pup implies Mc 4 x 
10~^ Mq yr~^, and is similar to the values suggested by the 
pulsation period and the K amplitude. 

The very high reddening value obtained for phase III 
supports our suggestion (Section 3.1 and 4.3) that the de- 
cline in the IR flux was due to obscuration of the Mira by 
dust. A much shorter phase of increased reddening also oc- 
cured in 1983/84. The possible causes of these events will 
be discussed later, however, the fact that the emission line 
ratios and the UV continuum do not indicate any increase in 
the reddening, implies that the site of the obscuration must 
be in the vicinity of the Mira leaving the hot component and 
the ionised region outside. 

After ~ JD 2 445 000 the K magnitude of RX Pup does 
not seem to be much affected by the hot component radi- 
ation, and can be used to estimate a the distance via the 
period luminosity relation (Feast et al. 1989; van Leeuwen 
et al. 1997): 

Mk = 0.88 - 3.47 log P. (3) 

Thus Mk ~ —8.7, although we note that this relation is 
not well calibrated for long period Miras. The average K 
magnitude corrected for the respective values of Eb-v did 
not change in epoch II-IV; < Ko >~ 2.64, and implies a 
distance of d = 1.8 ± 0.5kpc. 

In summary, the above analysis is consistent with con- 
stant reddening, Eb-v ^ 0.8 ± 0.2 mag, towards the hot 
component and the ionised nebular region, with an addi- 
tional variable reddening of Eb-v ~ to 2.5 mag towards 
the Mira. In the following analyses of the hot component 
and the ionised nebula values of Eb-v = 0.8 and d = 1.8 
kpc are adopted. 

3.3 Temperature and luminosity of the hot 
component 

Three methods are used to determine the temperature, Th, 
of the hot component. Wherever possible, Th is estimated 
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from the observed equivalent width of Hell A1640 emis- 
sion line from IlJE spectra, and by the Hell A4686/H/3 
or He 11 A4686/He I A5876 ratios from optical spectra. Since 
H/3 may be weakened by optical depth effects, we assume 
Hell A4686/He I A5876 is a more reliable indicator of the 
temperature than HenA4686/H/3 for the period 1983-1990, 
when the optical depth effects seem to be important. The 
third method applies the relation between Th and the high- 
est observed ionisation potential, Xmax, 

Th [lOOOK] = Xmax [eV], (4) 

proposed by Miirset & Nussbaumer (1994). The accuracy of 
this method is about 10 per cent provided that the highest 
ionisation stage is indeed observable. The highest observed 
ionisation stages as well as the derived temperatures at var- 
ious epochs are presented in Table 6. The highest ionisation 
stage observed in lUE spectra of RX Pup is that of Mg^'*, 
and in optical spectra are those of Fe"*"^ and O^^; they sug- 
gest temperatures in good agreement with estimates based 
on other methods. 

To derive the luminosity of hot component (s), Lh, the 
following methods are used: 

(a) The optical and IR photometry obtained during visual 
maximum (~ JD 2 442 500 - 4 000) suggests negligible con- 
tribution of the cool giant to the UBV magnitudes. While 
the JHKL magnitudes contain contributions from both the 
cool and the hot components. We assume that the IR mag- 
nitudes observed in epoch IV (see Section 3.2.4) represent 
those of the Mira, and that the contribution of the hot com- 
ponent to the IR continuum was negligible at that time. To 
estimate the JHKL magnitudes of the hot component in 
remaining epochs, we simply subtracted the contribution of 
the Mira from the observed values. The values for Lh in the 
last column of Table 6, correspond to the integrated (UV)- 
optical-IR luminosity. We combined the UBV magnitudes 
for JD 2 442 524 (Feast et al. 1977) with the JHKL data 
for JD 2 442 520, and the UBV magnitudes for JD 2 443 627 
(Klutz 1979), with the lUE spectra taken on JD 2 443 770 
and the average IR data from epoch I. 

(b) Our second method is based on the flux of the 
HellA1640 emission line measured on lUE spectra, or the 
He II A4686, H/3 and Ha fluxes derived from optical spectra. 
We assumed a blackbody spectrum for the hot component, 
and that Hi and Hell emission lines are produced by pho- 
toionisation followed by recombination (case B). The case B 
assumption is apparently not satisfied for the Hi emission 
lines in 1983-1990, and thus the corresponding values for Lh 
in Table 6 are those derived from Hell. Or, for JD 2 447 510, 
from Hof, whose flux should not differ from its case B value 
by more than 20 per cent (Netzer 1975). The HellA1640 
fluxes measured on low resolution lUE spectra before 1983 
have been corrected for the contribution of the Ol] A1641 
line estimated for the high resolution spectra. This contribu- 
tion was ~ 40 per cent of the total (He ll A1640 + O l] A1641) 
flux in 1980 September, ~ 30/25 per cent in 1981 June/1982 
March, and it becomes negligible after 1983 October. 

(c) Mikolajewska et al. (1997) found a relation between the 
hot component luminosity in symbiotic binary systems and 
the flux emitted in the Raman scattered O vi A6825 emission 
line. Schmid & Schild (1990b) report a weak emission at 
A ~ 6830A on 1 January 1990, with an equivalent width 



of ~ 11.7 A which they identify as Ovi A6825. Similarly, 
Harries & Howarth (1996) estimate EW(A6825) ~ 0.15 ± 
.05 A in May 1994. We converted these equivalent widths 
to fluxes estimating the continuum level from published Rc 
photometry (Munari et al. 1992; IS94). O vi A6825 was below 
our detection limit, ~ 10~^^ ergs~^ cm~^, on all our optical 
spectra, including those taken in November 1990. 

(d) In a few cases Lh was derived from UBV magnitudes 
observed during the nebular phase. The method was pro- 
posed by Miirset & Nussbaumer (1994), who also gave the 
bolometric corrections to the UBV magnitudes for a wide 
range of Th. For JD 2 447 978 and JD 2 448 215/23 published 
UBV magnitudes (Munari et al. 1992; IS94) were used, while 
for JD 2 448 328, the UBV magnitudes were calculated from 
our HST spectrum. 

(e) The HST spectrum was also used to calculate the in- 
tegrated (AA1770 - 5450A) flux. The flux emitted short- 
wards of 1200 A, i'Euv.h, can be estimated from the optical 
H/9, HelA5876, and HellA4686 emission lines (Kenyon et 
al. 1991; Eq.(8)). The sum of the integrated HST flux and 
i'Euv.h derived from the data in Table 3 for JD 2 448 352 
and JD 2 448 385 agrees very well with the values indicated 
by the method (b), and is by a factor of ~ 2 lower than the 
value derived from UBV magnitudes (method d). 

We estimate the luminosities in Table 6 have accuracies of a 
factor of 2, as indicated by the values obtained at the same 
epoch, but using different methods. In addition, there is a 
systematic error due to uncertainties in reddening and dis- 
tance estimates; this, however, should not affect any conclu- 
sions about the evolution of the hot component luminosity 
with time. 

The data in Table 6 demonstrate that the hot compo- 
nent had roughly constant luminosity in 1975-1989, although 
its temperature changes by a factor of 15. Between 1989 May 
and 1991 April, Lh dropped by a factor of ~ 20, while Th 
decreased by ^ 40 per cent. The evolution of RX Pup in 
1975-1991 thus parallels the evolution of a typical symbiotic 
nova (compare Mikolajewska & Kenyon 1992b; Kenyon et 
al. 1993; Miirset & Nussbaumer 1994). 

In 1995-96, the optical continuum is much stronger than 
in 1991, while the nebular emission is significantly reduced. 
There is no evidence for high ionisation lines. In particular, 
we do not see HellA4686 in a spectrum taken in 1996, and 
the upper limit for the equivalent width of the Raman scat- 
tered Ovi A6825 emission line, EW(A6825) <, 0.15A is com- 
parable to the value measured in 1994 by Harries & Howarth 
(1996), and requires Lh of order of 100 Lq. On the other 
hand, RX Pup was detected by ROSAT in April/ June 1993 
(Murset, Wolff & Jordan 1997), afihough the /C/i; spectrum 
from May 1993 does not show any evidence for high ion- 
isation features, nor does the optical spectrum from June 
1993 (IS94) which is very similar to our 1995/96 spectra. 
The lack of measurable Hell emission (for Th <; 10^ K and 
Lh ~ 100 Lq the observed flux in Hell 4686 should be of 
~ 10"^'' erg s~^ cm~^) may then be due to obscuration of 
the He ll formation region, which must be located very close 
to the hot component. 

The Hi Balmer emission lines are relatively strong, 
and their equivalent widths, EW(Ha) « 200 - 250A and 
EW(H/3) ?a 16 - 24A, measured on our 1995-96 spectra, 
respectively, are comparable to those observed in quies- 
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Table 6. Hot component parameters for RX Pup 



JD 24... Ion Th [K] [L©] Method 

a 
a 
b 
b 
b 
b 
b 
b 
b 
b 
b 
b 
b 
b 
b 
c 
d 
b 
d 
d 
e 
b 
b 
c 



cent recurrent novae T CrB and RS Oph (Dobrzycka et al. 
1996a; Anupama & Mikolajewska 1999). The He I emission 
lines are, however, very faint or absent. The colour tem- 
perature of the optical AA4500 — 5500A continuum is low, 
T - 7000 - 6000 KforEB-v = 0.8. This component has a 
luminosity of L ~ 6OO(d/1.8kpc)^L0. On the other hand, 
the optical emission lines indicate a higher temperature, 
T ~ 2.5 - 4 X 10" K, with a roughly comparable luminosity, 
L ~ 500 - 1000(d/1.8 kpc)^L0. We find here a striking sim- 
ilarity between RX Pup and the quiescent symbiotic recur- 
rent novae RS Oph and T CrB, as well as the hot component 
of CH Cyg during bright phases. In all these systems, the op- 
tical data indicate the presence of relatively cool F/A- (CH 
Cyg) or A/B-type source (RS Oph and T CrB), while the 
optical emission lines indicate a higher temperature source, 
with a roughly comparable luminosity (Mikolajewska et al. 
1988; Dobrzycka et al. 1996a; Anupama & Mikolajewska 
1999). We will return to this problem in Section 4. 

3.4 Nebular emission 

The rich emission-line spectrum that has been developing 
since late 1979, offers various possibilities for the determi- 
nation of the nebular parameters Te and ric- Published es- 
timates for the electron density indicate systematic differ- 
ences between the values derived from the UV intercombi- 
nation lines and those derived from the forbidden optical 
emission lines. In particular, K82 estimate Ue ~ 10^" cm~"^ 
from lUE observations in 1981 June, while W84 derive 
ne ^ lO'^ and ~ 10^ cm-^ from the [S ii] A6717/A6731 
and [Olll] A4363/(A4959+A5007) line intensity ratios, re- 
spectively, observed in 1983. Based on the optical spectrum 
taken in December 1989, FPC92 simulate the density gradi- 
ent inside the ionised region by adopting a two-layer model. 



in which different ionisation stages are present, and derive 
Tic ~ 4.5 X 10*^ cm~^ from high excitation [Olll] and [Fevi] 
emission line ratios, and Ue ~ 2.6 x 10® cm~^ from low ex- 
citation [N ll] and [S ll] line ratios, respectively. All these 
estimates assume a photoionised region with Te ~ 10* K. 

The significance of the density gradient in the ionised 
region is also indicated by the O ill] /[O ill] line ratio. Nuss- 
baumer & Storey (1979) demonstrated that simultaneous 
observations of the UV intercombination O ill] AA1661, 1667 
and the optical forbidden [O ill] AA4363, 4959, 5007 emission 
lines enable both Te and to be determined. Although 
the lUE and optical observations presented here are not 
simultaneous, we can combine the O ill] A1664 fluxes from 
Table 2 with published measurements for the optical [O ill] 
lines. In particular, the 1983 [O ill] fluxes from W84 com- 
bined with the mean of O ill] A1664 flux measured on the 
lUE spectra taken in 1982 and 1983, require a region with 
rie ~ 10^ - 10^° cm"^ and Te - 7500 K for ^b-v ~ 0.8 - 1. 
The Olll]/[Olll] intensity observed in the late 1980s (Ta- 
ble 2; van Winckel et al. 1993; FPC92) also require high 
densities and very low temperatures Te ^ 7000 K. Such a 
low temperature is, however, in conflict with the absence 
of strong lines arising from dielectronic recombination. A 
similar effect is seen in other symbiotic systems, including 
symbiotic novae (e.g. Schmid & Schild 1990a). To solve this 
problem, Schmid & Schild inferred that the forbidden lines 
come from a much more extended region (by a factor of 10* 
or so) than the O ill] emission, both having similar Te Jl 10* 
K. 

There is no simple method of determining Te in sym- 
biotic systems. Nussbaumer & Storey (1984) have proposed 
the use of dielectronic recombination lines. The idea rests 
on a comparison of a collisionally excited line with one pro- 
duced by dielectronic recombination. Nussbaumer & Storey 
also list promising recombination lines longwards of 912 A. 
The best candidates among these are Niv A1719 and Cll 
A1335. The observed flux ratio N v A1240/N iv A1719 grew 
from fa 1.5 in 1981 to 3 in 1984, and ?a 8 in 1986. The cor- 
responding temperatures are Te ~ 12 000 K in 1981, 13 000 
K in 1984, and 15 000 K in 1986, for ^b-v ~ 0.8. A very 
weak C11AI335 line is present on spectra taken in 1979 and 
1980. The observed flux ratios C ill] A1909/C 11 A1335 ~ 60 
in 1979, and ^ 100 in 1980 indicate Te ~ 12 000 - 13 000 
K. These temperatures argue very strongly in favour of a 
radiatively ionised UV emission line region. 

The fine ratios within the UV O iv] A 1401 and 
Nlll] A1750 multiplets are very sensitive to Ue (Nussbaumer 
& Storey 1979, 1982). The Nlll] A1754/A1749.7 and Nlll] 
A1752.2/A1749.7 ratios measured on lUE high resolution 
spectra indicate the electron density was decreasing from 
Ue 10^° - 10^^ in 1980 to ~ 5 x 10^ in 1982-83 and 
r-^ 10^ cm"^ in 1987. Similarly, we estimate Ue ~ 10^" in 1981 
and ~ 10^ cm-3 in 1987 from the O iv] A1407.4/A1404.8 ra- 
tio. These results are practically independent of Te. We also 
observe a decreasing trend in the Si ill] A1892/C ill] A1908 
and an increasing trend in the N ill] A1750/C ill] A1908 line 
ratios measured on the low resolution lUE spectra. These 
trends can be interpreted as a continuously decreasing Ue- 
Finally, the increase in the C iv doublet ratio in 1980- 
1988 (see Section 3.1.2) is consistent with decreasing optical 
depth (which is related to density) in the high-excitation 
permitted line region. 
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We can now determine whether or not the parameters 
derived from the intercombination lines are also representa- 
tive of the Hi/ Hel emission region. The ionisation poten- 
tials of 0+^ (54.9 eV) and He+ (54.4 eV) are very similar, 
so He I and O lll] should be emitted in the same physical re- 
gion. The fluxes in Hel A5876 and Olll] A1664 then require 
rie ~ 10^° cm"^ in 1983-84 and - 10^ cm"^ in 1988/90, 
and Te ~ 15 000 K in 1984-90. The dependence of emis- 
sivity in Hel A5876 on Te and We was taken from Proga 
et al. (1994) and Proga (1994), while for Olll] A1664 data 
from Nussbaumer & Storey (1981) were used. The volumes 
needed to account for the Olll] A1664 and Hel 5876 emis- 
sion fluxes are ~ 8 X 10"° and ~ 2 x 10"^ (d/l.Skpc)^ cm^ 
in 1983-84, and 1988-90, respectively, so the corresponding 
radii are i? ~ 2 and ~ 5 a.u. The resulting emission measures 
are a factor of ~ 2 higher than corresponding values derived 
from the H/3, and roughly equal those derived from the Ha 
emission fluxes. However, these estimates assume that the 
Hi lines follow case B. For a constant density gas, the opti- 
cal depth in Ha is thc ~ 10(7?/1 a.u.)/(ne/10^° cm"^) (Cox 
& Matthews 1969). We estimate th^ ~ 20 and ~ 5 in 1983- 
84 and 1988-90, respectively. In such conditions the H/3 flux 
could be reduced by ~ 50 — 70 per cent with respect to the 
case B value, while Ha should be somewhat (by ^ 20 per 
cent) enhanced. These results indicate that the UV inter- 
combination lines form in the same gas as do the H I and 
Hel lines. 

Our spectra taken in 1991 do not show any measur- 
able UV emission lines, while the optical Hi, Hel emission 
lines remain strong. The Hel line ratios are consistent with 
formation in a region with rie ~ 10® - 10* cm (Proga et 
al. 1994), while the Hi line ratios are consistent with case 
B for Es_v ~ 0.8 — 1. We estimate the emission measure 
nlV ~ 1.4 X lO^'* and ~ 1.2 x 10^^ (d/1.8 kpc)^ cm-^ for 
the Hel/Hl line emitting region (Te ~ 15 000 K), and the 
expected O lll] A1664 flux of ^ 7 x 10"^^ erg s"^ cm"^ which 
is below the detection limit on our lUE spectrum. 

We now examine the low density region in which 
the forbidden lines are produced. Based on an analysis 
of the forbidden lines in three symbiotic novae Schmid & 
Schild (1990a) found that physical conditions vary strongly 
throughout the nebulosities. In particular, they found a 
steep electron density gradient from the lowest to the high- 
est observed ionisation stages. Following their suggestion, 
we have grouped the diagnostic line ratios used according to 
the ionisational potential of the ions involved: 

(a) the transition region between the ionised and neutral 
gas, characterised by [Sll] and [Ol]; 

(b) the He^ zone and the outer edge of the He"*"^ region, 
described by [Nil], [Felll], [Arm], [O lll] and [Neiv]; and 

(c) the inner He"*"^ region, characterised by [Nev], [Mgv], 
[Fevi] and [Fevii]. 

The [S 11] A6716/A6731 ratio is consistent with We ~ 
few x 10" cm-^ and Te ~ 10* K (Canto et al. 1980) for 
the entire observing period. The electron temperature can 
be in principle constrained by [O l] (A6300+A6364)/A5579 
(Keenan et al. 1995), and the data in Table 3 are consistent 
with We ~ 2 x 10* cm"^ and Te ~ 10 000 K in 1991, and 
We ~ 5 x lO^cm"^ and a rather high Te <, 18 000 K in 
1995. However, [O l] A5579 can be severely affected by the 



atmospheric airglow (Keenan et al. 1995) so the [O l] ratio 
only constrains an upper limit for Te- 

The strong [Ar ill] A7135 and very weak or absent [Ar ill] 
A5192 lines (Table 3; FPC92) also indicate low Te. Our es- 
timate of /(A7135)//(A5192) Z 100 from the 1991 spectrum 
implies Te 15 000 K for Eb-v ~ 0.8 (Keenan, John- 
son & Kingston 1988). Unfortunately, our 1990 spectrum 
does not cover the ~ A 5200 A region. The absence of [Ar ill] 
A5192 in Table 2 of FPC92 suggests Te 20 000 K in 1989. 
The [Nil] A6584/A5755 ratio observed in 1989-91 indicates 
Ue ^ 1.5 X lO^cm-3 for Te <, 15 000 K (Kafatos & Lynch 
1980). Similarly, we estimate ^ 8 x 10*^, ^ 4 x 10*^, and ^ 
2 X 10® cm"^ from [Olll] (A5007-|-A4959)/A4363 observed in 
1983 (W84), 1989 (FPC92) and 1991 (Table 3), respectively 
(Nussbaumer & Storey 1981), and Ue ~ 5 — 9x 10® cm~^ from 
the [Neiv] A1602/A2423 ratio measured in 1986-87 (Nuss- 
baumer 1982). The [Felll] intensities are consistent with 
He ~ 10^ - 10® in 1990, and He ~ lO'^ cm"^ in 1991, re- 
spectively (Keenan et al. 1993). The volume needed to ac- 
count for L([Oiii] A5007) is V ~ 1.7 x lO"®, ~ 1.6 x lO"®, 
and ~ 6 X lO"^ (d/1.8 kpc)^ cm^ in 1983, 1988-90 and 1991, 
respectively, (Nussbaumer & Storey 1981), and the corre- 
sponding radii are R ~ 107/103/75 (d/1.8kpc)^/^ a.u. for a 
spherically symmetric nebula. 

There is also no evidence for high Te in the region where 
the emission lines with the highest ionisation potential orig- 
inate. [Mgv] AA2784, 2930 are very strong on several lUE 
spectra taken in 1986-89; A 2417 is, however, very weak or 
absent. Our estimate of J(A2784 -I- A2930)/7(A2417) <; 30 
is consistent with Te ^ 20 000 K for Ue ^ 10* cm"^, and 
FiB-v ~ 0.8 (Kafatos & Lynch 1980). Similarly, the absence 
of [Arv] A4626 on our optical spectrum taken in 1990, im- 
plies the observed ratio [Arv] 7(A6435 -h A7005)/7(A4626) ^ 
6, and Te <, 20 000 K for Ue <, 10* cm"^ and Eb-v ~ 0.8. 
The [Fevi] A4967/A4973 and [Fevii] A4967/A4989 ratios 
measured on our 1990 spectrum are then consistent with 
Ue ~ 10® -10^ cm~^ (Nussbaumer & Storey 1978; 1982). The 
volume needed to account for the luminosity in the [Fevil] 
AA6086, 5721 lines is then comparable to that required for 
[O III] lines if Fe/H is similar to that found in other symbiotic 
novae and M giants (e.g. Schmid & Schild 1990a). 

The physical conditions in the nebula are thus con- 
sistent with those expected from a photoionised gas. The 
bulk of permitted and intercombination lines are formed in 
a dense region of ~ a few a.u. and Ue ~ 10^ — 10^° cm~"^, 
with an average ~ 15 000 K. The forbidden lines come 
from more extended region of ~ 100 a.u. ([Olll], [Fevil]) 
or larger. We also find a steep density gradient covering 
more than 2 orders of magnitudes from the highest (Fe^®) 
to the lowest (S"*") observed ionisation stage in the forbidden 
line emitting region. Te also decreases with decreasing ion- 
isation. The conditions in the nebula were evolving during 
the period under discussion. The average Te was increasing 
with increasing ionisation level till c. 1989, then in 1991 it 
decreased in step with the ionisation degree. The opposite 
trend is observed in He; it was slowly decreasing till c. 1989. 
In both the permitted/intercombination line emitting region 
and in the forbidden line region the average Ue decreased by 
an order of magnitude. This trend is also manifested by a 
general increase in the forbidden emission line fluxes with 
respect to the permitted lines with similar ionisation poten- 
tials. 
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A comparison of the conditions derived from the Hi, 
He I and Olll] lines with those derived from the forbidden 
lines, [O ill] and [Fe vil] , suggests the density gradient ob- 
served in RX Pup is roughly of a form rie ~ r~^, which is 
consistent with the emission from an ionised wind emanat- 
ing from the system. A similar conclusion is reached from an 
analysis of the radio observations (HoUis et al. 1986; Seaquist 
& Taylor 1987, 1992). This shows the spectrum and angu- 
lar structure expected from partially optically thick ther- 
mal bremsstrahlung originating in a spherically symmetric 
steady stellar wind. The cm-mm spectrum observed in 1985- 
1988 is essentially linear from 1.49 to 394 GHz, and is well 
represented by a spectrum of the form 5*^ = 9.9 u"'^'^ mjy 
(Seaquist & Taylor 1992). The photospheric radius at 394 
GHz is given by r ~ 29(d/1.8kpc) a.u., and the turnover 
frequency ut 394 GHz implies ne <; 8 x lO'^ cm~^ inside 
that radius. 

For a spherically symmetric steady wind emanating pre- 
dominantly either from the cool giant or the hot companion, 
we have: 



where r is the radial distance from the centre of the 
star, the mean atomic weight of the gas and mn the 
mass of a hydrogen atom. Using the characteristic radii 
and densities derived above from the emission lines and 
radio mm data for 1985-88, we estimate A'l/v ~ 6 x 
10~^ [M© yr~^/100kms~^]. FW's analysis of emission line 
profiles indicate v ~ 140kms"\ and M ~ lO"'^ M© yr"\ It 
is unlikely that the cool giant wind has a velocity in excess 
of 100 kms~^. We thus assume the wind originates from the 
hot star. The evolution in the nebular conditions can then 
reflect the changes in that wind. Note that our estimate 
for M, although acceptable for a symbiotic nova (e.g. Kato 
1997), is very approximate. In particular, there is evidence 
for elongated structure in the radio and in the optical emis- 
sion lines (HoUis et al. 1986; Paresce 1990), which indicates 
the wind may be bipolar. The observed radio spectral in- 
dex of ~ +0.8 also departs from the canonical value of -1-0.6 
expected for a spherically symmetric wind. 

The change in the emission spectrum in 1991 can 
be understood in terms of decreased M and gradual re- 
combination in the ionised region following the decline in 
the hot component luminosity, which is no longer able 
to maintain the ionisation of the wind. As the recombi- 
nation time scale. Tree, is inversely proportional to Ue - 
Troc ~ lhr/1 month/ 15 yr for ~ lOVlOVlO* cm-^ re- 
spectively - the inner portions of the winds are first to re- 
combine. In fact, spectra taken in 1995 and 1996 show low 
ionisation emission lines of [O l] , [S ll] , [N ll] with almost the 
same (within a factor of ~ 2) intensity as observed in 1990 
and 1991, and very faint [O ill] lines; all these lines have 
intensities roughly consistent with their originating in a re- 
gion with He ~ Wcm"^ and Te 10 000 K. Such a sce- 
nario is also consistent with the observed changes in radio 
emission. This suggested an optically thick wind emanating 
from the system in the 1980s; while since c. 1991 it is better 
understood as optically thin emission (IS94, and references 
therein) arising from the low density remnant of the 1980-90 
high-ionisation phase. 

Only the Balmer Hi and Fell emission lines show in- 



creased intensity on our 1995/96 spectra. This increase coin- 
cides with the appearance of an A/F-type continuum source. 
The very steep Balmer decrement is consistent with its origin 
in a relatively dense region, Ue J5 10^" cm~'^ (Drake & Ul- 
rich 1980). Assuming Te 10* K, we estimate the optically 
thin emission measure n^V ~ 4 x 10^®(d/1.8kpc)^ cm~^ 
from the H/9 and Ha emission line fluxes. Again, we em- 
phasise here the similarity to the spectroscopic appearance 
of CH Cyg during its active phases. In particular, the lumi- 
nosity of the A/F-type continuum, as well as the conditions 
(density, emission measure) in the Hll Balmer emission re- 
gion at maximum of CH Cyg's activity (Mikolajewska et 
al. 1987, 1988) are comparable to those derived from the 
1995/96 spectra of RX Pup. 



4 DISCUSSION 

4.1 Geometrical configuration of RX Pup 

The large difference between the reddening towards the Mira 
and the reddening towards to the hot component and ionised 
region indicates the hot component and most of the ionised 
material must lie along a line of sight that does not include 
the Mira variable and its associated dust cocoon. 

Kenyon, Fernandez-Castro & Stencel (1988) deduced a 
characteristic dust radius of ~ 200 a.u. from IRAS colours, 
while Anandarao, Taylor & Pottasch (1988) derived a dust 
shell radius of ~ 77 a.u. (re-scaled for our adopted distance 
of 1.8 kpc) from analysis of the IRAS photometry and low 
resolution spectroscopy. Both values are an order of mag- 
nitude larger than the binary separation, a ~ 17 a.u., sug- 
gested by the IS94 model. 

Schmid & Schild (1990b) report weak emission at 
A 6825 A in 1990 January due to Raman scattering of the 
Ovi A 1032 resonance line. They also noted that the lack 
of detectable polarization in this line indicates the line of 
sight coincides roughly with the direction defined by the 
O"""^ zone and the Mira variable. RX Pup shows at the same 
epoch very strong [Fevil] emission lines, and the lack of 
strongly scattered O vi A 6825 emission is unusual. Practi- 
cally all symbiotic systems with strong [Fe vil] lines have also 
strong Ovi AA6825, 7082 lines (Kenyon 1986; Mikolajewska 
et al. 1997). The apparent weakness of the Raman scattered 
O VI line suggests that the scattering process in the envelope 
of RX Pup must be very inefficient, probably due to insuf- 
ficient numbers of the scatterers - neutral hydrogen atoms. 
Theoretical simulations of the Raman scattering in symbi- 
otic systems demonstrated that the process has the highest 
efficiency in systems having ionisation geometry with an X- 
parameter (as defined by Seaquist, Taylor & Button 1984, 
hereafter STB84), X ~ 1, while in systems with X ^ 5, 
the efficiency is much lower (Schmid 1996). The weak Ovi 
A 6825 would thus suggest that in RX Pup only a cone- 
shaped region shielded by dense material surrounding the 
Mira remains neutral, while the lack of measurable polar- 
ization indicates that in 1990 the 0+^ zone (and the hot 
component) was in front of this neutral region. Moreover, 
the lack of changes in the reddening towards the hot com- 
ponent and the nebular region (Section 3.2) shows that at 
least during the last 18 years both the hot component and 
the emission line region have remained outside the Mira's 
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dust shell. In particular, this rules out the possibility that 
the changes in the extinction towards the Mira component 
are orbitally related (we return to this problem in Section 
4.3). 

The conical shape of the interface between neutral and 
ionised regions in RX Pup can be produced either by illu- 
mination of the Mira's wind by the hot companion (STB84) 
or by interacting winds (Girard & Willson 1987). We have 
already demonstrated that the nebular emission in RX Pup 
originates from an ionised wind emanating from the system, 
and that the hot component is the most likely source of 
this wind. One can expect that this wind will interact with 
the Mira's wind. Assuming that the shock front coincides 
with the interface between the ionised and neutral regions, 
then X ^ 5 (see above) corresponds to a momentum ratio, 
mv = Mi,Vh/McVc 10. Adopting Mh ~ lO'^^MQyr'^ 
Vb ~ 140 kms"^ (Section 3.4), Mc ~ 4 x 10"® Mq yr'^ (Sec- 
tion 3.2.4), and assuming reasonable Vc ~ 20kms~^ (White- 
lock et al. 1994), we estimate mv ~ 18, in agreement with 
the suggested geometry. The distance of the shock from the 
Mira normalised to the binary separation will then be ~ 0.2. 
The latter result combined with the fact that the Mira com- 
ponent is never stripped of its dust envelope (Table 5) im- 
plies a binary separation that is a factor of ~ 5 or so larger 
than the dust formation radius. Assuming a typical dust for- 
mation radius of J; 5 x 7?c, and the Mira radius, _Rc ~ 2 — 3 
a.u. (e.g. Haniff, Scholz & Tuthill 1995), the required binary 
separation is a <; 50 a.u. 

With so large a binary separation, the Mira in RX Pup 
should have SiO masers, whereas the search for such emis- 
sion in 1987 (43 GHz; Allen et al. 1989), 1989 and 1994 (86 
GHz; Schwarz et al. 1995) resulted in an upper limit which 
was a factor of ~ 100 below fluxes in single Mira variables 
with similar pulsation periods. The popular explanations for 
the apparent absence of SiO masers (as well as OH and 
H2O) in most symbiotic systems involve radiation or tidal 
effects related to the hot compact companion (Schwarz et al. 
1995, and references therein). These are unlikely to apply to 
RX Pup, mostly because of the large binary separation, but 
also because a significant amount of the circumstellar dust is 
always present outside the SiO formation region (if the SiO 
molecules are destroyed by hot radiation the dust should be 
destroyed too). Moreover, at least the 1994 observation was 
made after the luminosity and temperature of the compan- 
ion significantly decreased. 

Allen et al. (1989) noted that if the cool giant were al- 
ways viewed through the ionised nebula, then the masers 
would be obscured at frequencies for which the nebula is 
optically thick. Moreover, Nyman & Olofsson (1986) found 
a correlation of the 87 GHz flux with the Mira pulsation 
phase. In fact, these two effects can explain the failure to de- 
tect SiO masers in RX Pup. First, our analysis suggests the 
Mira was always viewed through the ionised nebula. Thus 
in 1987, when the optically thick radio emission up to ^ 300 
GHz was observed (IS94) the 43 GHz maser would indeed be 
obscured. It seems that the radio spectrum remained opti- 
cally thick until c. 1991 (IS94), which indicates that in 1989 
the 87 GHz maser could be also have been obscured. Al- 
though the nebula was optically thin at 87 GHz at least in 
1994, both the 1989 and 1994 observations were made near 
the Mira minimum (0 ~ 0.42 and <j) ~ 0.48, respectively). 
The lack of detectable SiO masers in RX Pup is therefore 



consistent with the suggested geometrical configuration of 
this system. It also seems reasonable to repeat the search 
for such emission near to the Mira maximum, and when the 
system is in a low excitation stage with a flat radio spec- 
trum. 

4.2 Evolution of the hot component 

Fig. 11 shows the evolution of the hot component in the 
Hertzsprung-Russell diagram - for the luminosities and tem- 
peratures derived in Section 3.3 - along with the cooling 
curves for 0.64 Mq and 0.55 Mq white dwarfs, respectively. 
The hot component appears to maintain a constant bolo- 
metric luminosity, < Lh >= 15 000 ± IOOOLq, from 1975 
to 1986, though its temperature varies by a factor 15. In 
1988/89 it turned over in the HR diagram, since then our 
estimates suggest a decline in luminosity initially at roughly 
constant temperature. Between May 1989 and April 1991, 
Lh dropped by a factor of ~ 20, while Th decreased by ^ 40 
per cent. The evolution of RX Pup during its recent ac- 
tive phase thus parallels the evolution of a typical symbi- 
otic nova pointing to a thermonuclear nova eruption (com- 
pare Mikolajewska & Kenyon 1992b; Kenyon et al. 1993; 
Miirset & Nussbaumer 1994), and in particular it rules out 
the accretion model proposed by K85, with its recent devel- 
opments by IS94. Neither can the accretion model account 
for the emission line profiles or the evolution of the radio 
spectrum (see Section 3.4 for details). The evolution of Lh 
is also inconsistent with the thermonuclear-runaway model 
proposed by AW88, which assumes a very slow hydrogen- 
shell fiash with a plateau phase lasting decades. In particu- 
lar, the AW88 model does not predict any significant changes 
in Lh , as all photometric and spectral changes have been at- 
tributed to changes in the mass loss from the white dwarf. 

The estimated maximum bolometric luminosity, Lh ~ 
15OOOL0 falls within the range of plateau luminosities 
for other symbiotic novae (Mikolajewska & Kenyon 1992b; 
Miirset & Nussbaumer 1994). Applying the relationship be- 
tween plateau luminosity and white dwarf mass for accret- 
ing cold white dwarfs from Iben & Tutukov (1996, their 
equation 8), we estimate the mass of the white dwarf in 
RX Pup, mwD ~ 0.59Mq, while the Paczyriski-Uus rela- 
tion (Paczynski 1970; Uus 1970) gives mwD ~ 0.77Mq. The 
luminosity-mass relation during the plateau phase, however, 
depends significantly on the thermal history of the white 
dwarf: generally hot white dwarfs must have larger masses 
than the cool ones to reach the same luminosity during a 
hydrogen- burning phase (Iben & Tutukov 1996). Our esti- 
mates for the white dwarf mass of RX Pup should thus be 
considered as lower limits. 

RX Pup may be the first recurrent nova detected in 
a D-type symbiotic system. Published observations suggest 
RX Pup was visually bright in 1894, and its spectrum re- 
sembled that of rj Car, with strong H I Balmer emission lines 
(Pickering 1897). Three years later a high excitation nebu- 
lar spectrum emerged, with strong emission lines of Hell, 
[O III] , [Ne III] , in addition to the H I Balmer series (Cannon 
1916). After 1900, RX Pup entered a deep optical minimum 
(mpg = 14.1 mag in 1904-5; Pickering 1914) resembling that 
in 1991. The star recovered from the minimum in 1909, and 
from then up until 1924 it showed only small fluctuations 
around mpg ~ 12 mag (Yamamoto 1924). The high excita- 
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Figure 11. Evolution of tlie liot component in tlie HR diagram 
in 1975-1996. Dots indicate the 1975-1994 outburst cycle. The 
shaded square represents the results derived from H I emission 
lines, while the open circle locates the A/F-type component ob- 
served in 1995-96. The solid curves are the evolutionary tracks 
from Schonberner (1989), the dashed curve is a line of constant 
radius. 



tion spectrum was again observed in 1941 (Swings & Struve 
1941), while in 1949-51 the excitation level was moderate 
and the optical continuum was very weak (Henize 1976). 
In the late 1960s RX Pup remained very faint optically 
(mvis = 12.5, mpg = 13.5) and the high excitation features 
were absent (Sanduleak & Stephenson 1973). In 1969/70 the 
V magnitude brightened to 9.4 mag (Eggen 1973), and the 
recent outburst began. The historical records of RX Pup 
thus suggest that around 1894 it underwent another nova- 
like eruption. Unfortunately, the lack of observations in the 
1930s makes it impossible for us to ascertain whether the 
presence of a high-excitation spectrum in 1941 was also due 
to the eruption. 

The recurrence time for RX Pup, Tree ^ 80(40) yr, is 
comparable to that of two symbiotic recurrent novae, T CrB 
(80 yr) and RS Oph (22 yr), and it requires an accretion rate 
of the order of ~ lO"''M0yr"^ (Prialnik & Kovetz 1995). 
According to Bondi-Hoyle theory even in a wide binary sys- 
tem a white dwarf can accrete as much as 1 per cent of 
the giant's stellar wind (Livio 1988, and references therein), 
which in the case of RX Pup is just what is needed. The 
observed quiescent luminosities of the compact components 
of Mira AB (~ 2Lq; Karovska et al. 1997)0 and R Aqr 
(a few Lq; Michalitsianos, Kafatos & Hobbs 1980) spatially 
resolved wind-accreting binary systems agree with the val- 
ues predicted by the accretion theory. The component sep- 
aration of Mira AB, a ~ 70 a.u. (Karovska et al. 1997) is 



^ Their original values are scaled to the photometric parallax, 
118 pc (e.g. Van Leeuwen et al. 1997), which is in good agreement 
with the Hipparcos trigonometric parallax. 



comparable to that expected for RX Pup, while the quies- 
cent Lh in RX Pup is a factor of ~ 100 higher than that of 
Mira B in good agreement with more or less the same order 
of magnitude difference in the mass-loss rates of the Mira in 
RX Pup (Afc ~ 4 X 10"^ Mq yr"^ Section 3.2.4) and that of 
Mira A (Mc ~ 10"''MQyr"^; Reimers & Cassatella 1985). 

The observed maximum luminosity as well as the very 
slow evolution in the H-R diagram indicate a white dwarf 
mass significantly less than the Mh ~ 1.4 Mq usually ex- 
pected for recurrent novae (e.g. Webbink et al. 1987). As the 
estimated Lh is based largely on UV emission line fluxes, 
it will be underestimated if the nebula is not ionisation 
bounded. However, any underestimate is unlikely to be sufii- 
cient to increase the white dwarf mass above ~ 1 Mq . More- 
over, recent model calculations show that while recurrent no- 
vae require high mass-accretion rates, M ^ 10~* Mq yr~^, 
the constraint on the white dwarf mass is less severe (Prial- 
nik & Kovetz 1995). Short recurrence times can be obtained 
for white dwarf masses as low as ~ IMq. Consequently, the 
decline times of recurrent nova models are not necessarily 
short. Given all the uncertainties involved in the theoreti- 
cal modelling of symbiotic and recurrent novae (Sion 1997, 
and references therein), it seems plausible to suggest that 
RX Pup is a recurrent nova. 

The quiescent hot component of RX Pup shows strik- 
ing similarities with its counterparts in RS Oph and T CrB. 
In particular, the ROSAT spectrum of RX Pup observed 
in April- June 1993 (Miirset et al. 1997) is very similar to 
the ROSAT spectrum of RS Oph (Orio 1993); both show 
P type spectra (according to the classification of Miirset et 
al.) compatible with emission from an optically thin plasma 
with a temperature of a few 10® K. Miirset et al. suggest 
that such emission is due to hot shocked gas, and propose 
that the shock arises in the region where the strong and fast 
wind from the hot component crashes into the cool giant or 
its wind. Such an interpretation is, however, entirely incom- 
patible with the observed properties of both RX Pup and 
RS Oph at the time of the ROSAT observations. The collid- 
ing wind model requires very Lh ~ 10 OOOLq for RX Pup, 
and ~ 2 OOOLq for RS Oph, which is much higher that the 
observed Lh ^ IOOOLq in RX Pup, and ~ 100 - 600Lq 
in RS Oph (Dobrzycka et al. 1996a). Moreover, the optical 
spectrum of RS Oph at that time was very similar to the 
quiescent (1993-96) spectrum of RX Pup, showing only weak 
low-ionisation and low-excitation emission lines (Dobrzycka 
et al. 1996a). 

It seems thus reasonable to assume (as did Orio for 
RS Oph) that the X-ray emission arises in the accretion flow. 
Assuming that the ROSAT flux represents a significant frac- 
tion of the total flux in the 0.2-4 KeV range, a mass accretion 
rate Mace « 10"''(Mwd/0.7Mq)"^-**Mq yr"^ (Patterson & 
Raymond 1985) can be estimated. However, this should be 
considered a lower limit to Mace, because the X-ray luminos- 
ity for RX Pup given by Miirset et al. was derived assuming 
a much lower reddening, E_b-v = 0.55, than that derived 
above. In addition, local material might absorb a substan- 
tial fraction of the X-rays produced close to the white dwarf 
surface. 

The quiescent optical spectrum of RX Pup, in particular 
the presence of a variable relatively cool F-type continuum 
with strong H I Balmer emission lines, also resembles the 
optical spectra of RS Oph and T CrB. These stars show an 
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additional variable component in the optical-UV and some- 
times also at IR wavelengths. For example, in RS Oph the IR 
colours indicate the presence of an additional warm, ^ 7000 
K source (Evans et al. 1988), while Dobrzycka et al. (1996a) 
found an A-B type shell source with L ^ 100 — 600Lq, ac- 
companied by strong Hi and moderate Hel emission lines. 
Similarly, Selvelli, Cassatella & Gilmozzi (1992; also Bel- 
czynski & Mikolajewska 1998; Anupama & Mikolajewska 
1999) reported variable UV/optical continuum with L ~ 
40 - IOOLq in T CrB. Although both in the symbiotic re- 
current novae, T CrB and RS Oph, and in RX Pup, the aver- 
age luminosity of the B/A/F-type shell source is consistent 
with the accretion rate, M ^ lO~*M0yr~^, required by 
the theoretical models, the effective temperatures places the 
hot components far from the standard massive white dwarf 
tracks in the HR diagram. Simultaneously, the X-ray data 
suggest 1-2 orders of magnitude lower accretion rates. All 
three systems show similar quiescent behaviour: their hot 
components have highly variable luminosity and occasion- 
ally display blue-shifted absorption features (see also Do- 
brzycka et al. 1996a; Anupama & Mikolajewska 1999). 

It is also interesting that the shell spectrum appears 
only during the late decline from the nova outburst. In par- 
ticular, the optical/visual light curves from the outbursts of 
T CrB, RS Oph and RX Pup show more or less pronounced 
minima followed by a standstill (or secondary maximum) as- 
sociated with the appearance of the shell spectrum. We sug- 
gest the minima are due to a decline in the hot component's 
luminosity after it passes the turnover in the HR diagram. 
The strong hot component wind during the plateau phase 
prevents accretion onto the hot component. Following the 
decline in luminosity the wind also ceases, and accretion of 
the material from the cool giant can be restored. We believe 
that the shell-type features and variable "false atmosphere" 
observed at quiescence together with the H I emission line 
spectrum originate from the accretion flow. The observed 
variability could be due to fluctuations in the mass-loss rate 
from the cool giant, the shell spectrum becomes stronger 
as Mc increases and as a result Mace increases. Finally, we 
recall the behaviour of CH Cyg, a highly variable symbi- 
otic system containing an accreting white dwarf. The hot 
component luminosity in this system varies by a factor of 
10" (Mikolajewska 1994), between ~ 0.1 and ~ 300 L©. The 
brightening of the hot component is associated with the pres- 
ence of a flickering, which suggests the system is accretion- 
powered. The optical/UV spectrum of CH Cyg during bright 
phases resembles the quiescent spectra of the symbiotic re- 
current novae (see also Section 3.3), while the detection of 
hard X-ray emission (Miirset et al. 1997; Ezuka, Ishida & 
Makino 1998) indicates that the accreting component is a 
white dwarf. 



4.3 Obscuration of the Mira component 

The light curves of symbiotic Miras were discussed by 
Whitelock (1987, 1988). RX Pup is typical of these in that 
it shows significant long term variations (refereed to as "ob- 
scuration events" by Whitelock) in addition to the Mira pul- 
sation. Fig. 12 shows the J and L light curves after removal 
of the Mira pulsations. This is done by subtracting the best 
fitting 578 day sine-curve from the data as displayed in Table 
4 and Fig. 4. Only J and L are illustrated as they represent 
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Figure 12. J and L light curves after the removal of the Mira 
pulsation. 



the extremes of behaviour. During the first approximately 
2000 days low amplitude 578 day variations can be seen in 
anti-phase with the pulsation cycle; these are artifacts of 
the pulsation removal process introduced by the lower am- 
phtude of the Mira at that stage (see Section 3.2.4). The A J 
curve changes over a full range of about 2.5 mag. While the 
AL light curve shows a much lower amplitude, ~ 0.8 mag, 
trend, with a minimum at the same time as the AJ curve. 

Although the light curves in Fig. 12 seem qualitatively 
similar to the visual light curve (Fig. 4), the optical min- 
imum around JD 2 448 300 can be accounted for by the 
evolution of the hot component, in particular, the general 
decline in the emission line spectrum (the change in emis- 
sion line fluxes alone can account for a ~ 1/1.5 mag drop 
in V/B magnitudes, respectively; Section 3.1) following the 
decrease in its luminosity. In fact, the visual magnitudes are 
generally not correlated with the IR magnitude, although 
the hot component may contribute to the IR before ~ JD 
2 444 000 (epoch I/optical maximum; see Section 3.2.4). In 
particular, the visual light curve was practically flat dur- 
ing JD 2 445 000-8 200, while the IR light curves showed 
pronounced changes (Fig. 12). Similarly, the decrease in IR 
flux after ~ JD 2 449 500 is not accompanied by any related 
changes in the visual. This, and the fact that the changes 
in the reddening towards the Mira are not correlated with 
similar changes in the reddening towards the hot component 
and emission line formation region(s), strongly suggest the 
obscuration affects only the Mira. 

The physical nature of this phenomenon remains a mys- 
tery. In the specific case of the symbiotic star R Aqr the long 
term trends were explained as the effect of an orbitally re- 
lated eclipse of the Mira by pre-existing dust (Willson et 
al. 1981; Whitelock et al. 1983b). However, more recently 
Whitelock et al. (1997) have explained large amplitude er- 
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ratic variations in the light curves of carbon-rich Miras as 
the consequence of the ejection of puffs of dust by the star. 
The wavelength dependence of these two effects: an eclipse 
by pre-existing dust and the ejection of new dust, are very 
similar and it is difficult to determine which is occurring 
in RX Pup. From Fig. 12 it is clear that any eclipsing dust 
cloud would have to be very patchy, much more so than that 
responsible for the eclipse in R Aqr. 

Fig. 12 suggests a modulation with a time scale of 
~ 3000 days, while the reddening towards the hot compo- 
nent was practically constant over at least 13 years. It is hard 
to imagine a geometry of the binary system in which the cool 
component is eclipsed at least twice while the hot compo- 
nent is not eclipsed at all. Thus, for RX Pup, it seem unlikely 
that the obscuration is orbitally related. Now that data are 
available over a longer time base it seems that the obscura- 
tion events occur in many symbiotic Miras with too great 
a frequency to be generally associated with binary eclipses 
(Whitelock 1998). 

Although the most pronounced obscuration event coin- 
cides with the hot component passing the turnover in the 
HR diagram the obscuration is apparently not associated 
with this component. In particular, the dust does not form 
in the material ejected from the hot component. The rel- 
ative intensities of Siiv A 1394 -I- Oiv] A1403/Niv A 1487 
and Olll] A 1664/Nlll A 1750 should not vary with changes 
in Th, ric and ionisation. The lUE data in Table 2 show that 
in RX Pup these ratios are roughly constant until at least 
JD 2 447 667 which means that at that time O and Si are 
not being depleted due to condensation into dust grains (sili- 
cates). Similarly, there is not indication for carbon depletion 
in the ejecta. Moreover, the absence of strong O I during the 
late outburst phase and the steady decline in Hen emis- 
sion suggest the hot component ejecta lack a neutral zone. 
Similarly, Schmid & Schild (1990a) do not find evidence for 
significant depletion of Si and Fe in the ionised ejecta of 
V1016 Cyg and HM Sge, while Munari & Whitelock (1989) 
did not find any major changes of emission line strengths 
associated with the dust obscuration episode in HM Sge. 

The site of the dust responsible for the symbiotic Mira 
obscuration must be located in its circumstellar envelope. 
In a considerable number of Mira variables, the visual light 
curves exhibit complex variations from cycle to cycle as well 
as on long time scales. Although the long term modulations 
are more evident in C-rich Miras (e.g. R For, RU Vir), they 
are also present in the light curves of longer (^ 300"^) period 
0-rich Miras (S CMi, R Hya; Mattel 1997). Longer period 
Miras generally have more dust and it seems likely that the 
long term modulation of their light curves is related to the 
behaviour of this dust. Because they have more dust they 
are visually fainter and therefore less well studied at short 
wavelengths. It is therefore possible that these obscuration 
events have been easily recognised in symbiotic Miras simply 
because their light curves have been more extensively stud- 
ied. Whitelock's (1998) discussion suggests this is not the 
case, but more extensive studies are necessary to confirm 
this. 

It is not at all clear how, or if, the obscuration events 
affect the rest of the symbiotic system. It is possible that 
an enhancement in the Mira wind produces a coordinated 
increase in the accretion luminosity of the hot component 
as suggested by the rise in the optical continuum and H I 



emission lines after ~ JD 2 450 000 (Section 3.3). A similar 
relation between the obscuration of the Mira (and related en- 
hancement in its wind) and an increase in the hot component 
luminosity is probably observed in R Aqr (Mikolajewska & 
Kenyon 1992b). 



5 CONCLUSIONS 

RX Pup is a long-period interacting binary system consist- 
ing of a Mira variable pulsating with P ~ 578 days, and a 
white dwarf companion. The binary separation could be as 
large as a <; 50 a.u. (corresponding to Porb 200 yr) as sug- 
gested by the permanent presence of a dust shell around the 
Mira component. In particular, the Mira is never stripped of 
its dust envelope contrary to the model proposed by IS94, 
and even during bright (not obscured) IR phases the star 
resembles the high-mass loss galactic Miras with thick dust 
shells. 

The analysis of multifrequency observations shows that 
most, if not all, photometric and spectroscopic activity of 
RX Pup in the UV, optical and radio range is due to activ- 
ity of the hot component, while the Mira variable and its 
circumstellar environment is responsible for practically all 
changes in the IR range. In particular, the evolution of the 
hot component in the HR diagram as well as evolution of 
the nebular emission in 1970-1993 is consistent with a sym- 
biotic nova eruption, with the luminosity plateau reached in 
1972/75 and a turnover in 1988/89. The hot component con- 
tracted in radius at roughly constant luminosity from c. 1972 
to 1986; during this phase it was the source of a strong stellar 
wind and therefore could not accrete any further material. 
By 1991 the luminosity of the nova remnant had decreased to 
a few per cent of the maximum (plateau) luminosity, and the 
hot wind had practically ceased. By 1995 the hot component 
start to accrete material from the Mira wind, as indicated 
by a general increase of the optical continuum and Balmer 
H I emission. The quiescent optical spectrum of RX Pup re- 
sembles the quiescent spectra of symbiotic recurrent novae, 
while the hot component luminosity is consistent with vari- 
able wind-accretion at a high rate, Mace ~ 10~^MQyr~^ 
(f5i 1 per cent of Mc). RX Pup may be a recurrent nova; 
there is some evidence that a previous eruption occurred 
around 1894. The proposed binary model for RX Pup could 
be verified if fiickering of the accretion disk were detected at 
quiescence. Such fiickering has been observed in the symbi- 
otic recurrent novae RS Oph and T CrB during quiescence, 
and in the accretion-powered symbiotic systems CH Cyg and 
MWC 560 (e.g. Dobrzycka, Kenyon & Milone 1996b). 

Large changes were found in the reddening towards the 
Mira, these were accompanied by fading of the near IR flux. 
However, the reddening towards the hot component and 
emission line regions remained practically constant and was 
generally less than that towards the Mira. These changes do 
not seem related to the orbital configuration nor with the 
hot component activity. The IR dust obscuration episodes 
are best explained as intrinsic changes in the circumstellar 
environment of Mira possibly due to variable mass-loss rate. 
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